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Tests at Allied Chemical show 
Nitrogen Tetroxide ready for in- 
stant use even after 9 years of 
storage. Combustion tests with 
important liquid fuels and Nitro- 
gen Tetroxide indicated up to 
99% theoretical combustion effi- 
ciency. Additional advantages 
offered by NoOg are listed at right. 
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page Product Bulletin, a bro- 
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dling of Nitrogen Tetroxide”, 
and “Nitrogen Tetroxide as an 
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The Dynamics of Conducting Gases 
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flight, the use of conducting gases, namely plasmas, 

has become an important central problem. Among 
the various aspects of interest to engineers and scientists 
are the properties and the dynamics of plasmas as well as 
the interaction between the fluid dynamic and the 
electromagnetic fields. 
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Largest closed di die The forging of these 657 lb. bulkheads presented 


unusual problems—problems of die design, of forg- 
ing, and of handling. The successful solution of these 


forging ever made problems saved a costly assembly of many smaller 


: mg parts—costly both in terms of dollars and pounds. _ 


for 


This is a typical Wyman-Gordon accomplishment 
—the result of over seventy-five years of forging 
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Over 1000 Wyman-Gordon forgings 


pace advanced jet designs 


In saluting the great achievement of America’s 
commercial jet leadership, Wyman-Gordon takes 
great pride in having participated so extensively 
in the development of the aircraft industry from its 
very beginning. 


In World War I, the use of forgings was confined to 
engine parts and Wyman-Gordon made most of 


In the highly stressed parts of modern aircraft, there is no 


@ na 7" substitute for forgings and in complicated forgings of diffi- 
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cult alloys, there is no substitute for Wyman-Gordon quality 
and experience. We are prepared to serve you at the design, 
engineering and purchasing stages of your developments. 
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them. Today, from the huge main bulkhead 
forging at the left to the engine parts below, Wyman- 
Gordon forgings are in the wings, fuselage, landing 
gear and engines of the four newest commercial jet 
aircraft. These and a great variety of other vital 
components are forged from many types of ma- 
terials: aluminum, magnesium, titanium, steel, as 
well as many new uncommon alloys. 
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Laboratories 


in Canoga Park 


environment dedicated to 
technological research and development 


The new Ramo-Wooldridge Laboratories in Canoga 
Park, California, will provide an excellent environment 
for scientists and engineers engaged in technological 
research and development. Because of the high degree 
of scientific and engineering effort involved in Ramo- 
Wooldridge programs, technically trained people are 
assigned a more dominant role in the management of 
the organization than is customary. 

The ninety-acre landscaped site, with modern build- 
ings grouped’around a central mall, contributes to the 


academic environment necessary for creative work. The 
new Laboratories will be the West Coast headquarters 
of Thompson Ramo Wooldridge Inc. as well as house 
the Ramo-Wooldridge division of TRW. 

The Ramo-Wooldridge Laboratories are engaged in 
the broad fields of electronic systems technology, com- 
puters, and data processing. Outstanding opportunities 
exist for scientists and engineers. 

For specific information on current openings write 
to Mr. D. L. Pyke. a 
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"PIN skin missile tanks stiffened by internal pressure still 

may fail by local buckling if the loads are high enough or 
the pressure is not high enough. (A cylinder is considered to 
have failed by buckling when it is in such a state that the 
cylinder cannot sustain additional compressive stresses.) 
Unpressurized thin walled cylinders definitely may be de- 
signed by buckling. Sometimes such structures are missile 
structures. In the latter instance reliable data are mandatory 
in design. Much attention has been given in the literature to 
the problems of the buckling of unstiffened circular cylindrical 
shells under various loading conditions (1-52).1_ However, 
to date, reliable and simple design curves still are not avail- 
able. Therefore, a careful objective review of the known 
theories and test data has been undertaken. The purpose of 
the present paper is to present the results of the investigation 
in forms which are considered reliable and simple. Also, 
areas are pointed out in which there is a need for further test 


data. 


The greatest concentration of effort has been for this load- 
ing condition. The reasons for this emphasis, paradoxically 
enough, are the relative simplicity of testing for this type of 
loading and the accompanying need for explanation of the 
wide variance of the test results. Actually before the de- 
velopment of missile structures, the primary loading on air- 
plane monocoques was bending. The literature contains 
ample discussion of lack of agreement of theory and test for 
cylinders in compression. The problem is nonlinear with the 
strength of the structure strongly dependent on eccentricities 
of the structure. The problem therefore must be treated on 
an empirical and statistical basis. 

The final curves recommended for design must be of a high 
degree of reliability, yet not overly conservative. It has 
been argued that certain test data should not be used to es- 
tablish design curves. For example, Donnell’s test speci- 
mens were manufactured from shim stock materials and were 
of small diameter. However, other investigators also used 
small gages. Yet Donnell’s tubes were the smallest in di- 
ameter. Furthermore, Newmark tested large thick fabricated 

Presented at the ARS 13th Annual Meeting, New York, N. Y., 
Nov. 17-21, 1958. 

1 Numbers in parentheses indicate References at end of paper. 
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cylinders. His test data fall in about the same region as Don- 
nell’s or even lower. The edge condition of Newmark’s 
cylinders apparently was not a fixed edge condition. The 
conclusion that may be drawn here is that size or gage of a 
specimen in itself is not indicative of the reliability of the 
test data. One can build a bad large cylinder or perform a 
bad test on a large cylinder just as easily as on a small one. 
The question to be answered then is what tests should be con- 
sidered in arriving at design curves. Since, in actual prac- 
tice, conditions may not be as ideal as in the test laboratory, 
one should use all representative test data. Machined 
cylinders should not be considered because many actual eylin- 
ders are made up of spliced sheets. Cylinders tested to date 
had longitudinal splices. In the future it may be of advan- 
tage to test large cylinders, such as occur in missile structures. 

Another matter influencing the buckling of axially com- 
pressed cylinders is the uniformity of the axial state of stress 
around the circumference of the cylinder. The closer this 
state of stress is to a uniform state of stress, the higher the 
buckling load. It is of interest to note that the higher values 
of buckling stresses were obtained by investigators who used 
a central means of load application, i.e., a ball and socket set- 
up or a central strut. On the other hand, the test data ob- 
tained by use of plate-distributed loading are the lowest 
values in the test scatter. 

According to (22), a cylinder in axial compression may be 
considered long if a certain parameter Z is greater than 80. 
Thus for R/t > 100 and L/R > 1, the buckling stress would 
be independent of L/R. However, a correct interpretation of 
the test data reveals this conclusion is not correct. Where 
large deflection effects are important, eccentricities are 
important because the structure is more apt to jump from 
one state to another. The longer the cylinder is the greater 
the chance for out-of-roundness. Consequently, the param- 
eter L/R must be considered throughout the entire range 
of values of R/t. This effect is illustrated indirectly in Fig. 
la (20). Notice that the slope of the best fit curve is less 
than the 45-deg slope of the theoretical lines. This fact 
means that the buckling stress decreases with increasing 
length. More test data are needed for R/t larger than 2000 to 
confirm this trend. Some test data taken from (22) are 
shown in Table 1. All cylinders are so-called long cylinders 
(Z > 80). Notice the drop in allowable stress with increas- 
ing length. 
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~ a Fig. la Buckling test data and corresponding empirical curve for clamped axially compressed circular cylinders 


Table 1 Effect of length on the buckling of cylinders in 
axial compression, R = 8.75 
Specimen 
no. L R/t Ccr 
25(1) 9.5 2734 1739 
29(1) 9.5 2734 1648 
4(1) 21.5 2734 1176 
4(2) 21.5 2734 1148 
7(1) 21.5 2734 1097 
7(2) 21.5 2734 1119 
58(1) 21.5 2734 1250 
67(1) 21.5 2734 1222 


Based on the previous discussions, the curves presented in 
Fig. lb are recommended for design. Most of the test data 
are for cylinders with clamped edges. The curves should be 
applicable for buckling of simply supported cylinders except 
perhaps in the upper left-hand region. The allowable 
buckling stresses for long cylinders with large R/t are found to 
be less than those previously recommended for design. The 
curves shown in Fig. 1b are believed to be of high reliability— 
at least 90 per cent reliable. However, they are not based on 
a comprehensive statistical analysis because there are not 
enough data available to cover adequately the entire range of 


Pure Bending 


The test data for pure bending buckling of cylinders are less 
abundant than those for pure compression for reasons ex- 
plained previously. The maximum bending buckling stress 
M/xR*t exceeds the compressive buckling stress of a given 
cylinder. The theoretical reason for the increase is the fact 
that the compressive stress is not constant. Superimposed on 
this basic phenomenon is the effect of imperfections which 
again causes a wide variance in the test data. One might 
expect the test scatter to be less than in pure compression 
tests since the region of maximum compression is confined to a 
small portion of the cylinder cross section, thereby lessening 
the chance of imperfections precipitating buckling. However, 
there is another phenomenon connected with the bending of a 
cylinder which does not occur in compression tests. This 
effect is flattening or ovulation which is tantamount to eccen- 
tricity. The weakening effect of length discussed for the 
compression allowables also should be present for bending 
buckling for the reasons given previously. 

A suggested overall factor of 1.3 has been given for the 
ratio of the pure bending buckling stress to the pure compres-_ 
sion buckling stress. According to the results of the limited 
test data, this factor may be all right for certain ranges of 
values of R/t and L/R. However, for large R/t the factor 
may actually approach unity. More tests are needed in this 
region, as well as for long cylinders with smaller values of R/t. 
Some values of pure bending buckling allowables may be 
found from Fig. 5 by setting the pressure equal to zero. 


Pure Torsion 


The amount of test data for this case is plentiful. The — 
status of this loading condition is better than that for pure 
compression or pure bending as far as agreement between 
theory and test is concerned. The buckling behavior of the 
structure is less nonlinear and less sensitive to imperfections, 
although these effects are still present. The curves proposed 
for design are shown in Fig. 2. Notice the similarity between 
the shapes of these curves and the ones for pure compression 
shown in Fig. 1. The curves in Fig. 2 are of a high degree of 
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Fig. lb Compression buckling stress for unpressurized cylinders 
with clamped edges 
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Fig. 2. Torsion buckling stress for unpressurized cylinders with 
simply supported edges 


reliability. Strictly, these curves are supposed to be valid for 
cylinders with simply supported edges. However, they may 
be used for cylinders with clamped edges, except perhaps in 
their upper left-hand regions. Actually the test data are 
based on specimens whose ends were clamped, or close to 
being clamped. 
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External (Negative) Hydrostatic Pressure 


There are enough test data for this loading condition to 
draw design curves. Effects of eccentricities and nonlinear 
action are confined to so-called short cylinders (having low 
values of Z). The high reliability design curves are given in 
Fig. 3. These curves are for simply supported cylinders. 
For cylinders with fully clamped edges, it is suggested that the 
curves be used with the effective length of the cylinder reduced 
to 0.83 times the actual length. The curves in Fig. 3 may also 
be used for cylinders subjected to radial pressure, except per- 


haps in their upper left-hand regions. eae: (er 
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Fig. 3 Hydrostatic buckling pressure for cylinders with simply 
supported edges 


Fig. 4a Compressive buckling stress as a function of pressure 
hoop strain, L/R = 1 
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Each of Figs. 1 to 3 contains a reduced modulus F£, needed to 
account for inelastic behavior of the material when the 
buckling stress exceeds the proportional limit. Apparently 
Bijlaard was one of the first to derive these quantities. 
Krivetsky (5) plotted these results in certain diagrams to 
facilitate their use. More recently (9, 17) Gerard has de- 
veloped expressions for plasticity reduction. Lee and Ades 
investigated the reduced modulus for plastic torsion buckling 
Needham has suggested a value of »/EtE for the case of axi- 
ally compressed cylinders (27). This factor is the so-called 
von Karman reduced modulus. 

Donnell and Wan have related the onset of plasticity in an 
axially compressed cylinder to the severity of the imperfec- 
tions existing in the cylinder. For purposes of design, it is 
believed simpler to present the design curves as shown here- 
including the effects of imperfections only indirectly and using 
a reduced modulus to account for inelastic behavior. 


Combined Compression and Internal Pressure 


This loading condition is important in missile design. The- 
oretically, one would expect that if the pressure is made high 
enough the buckling stress could be made equal to the classi- 
cal value of 0.605 Et/R, since the internal pressure smooths 
the surface of the cylinder, thereby canceling the effects of 
imperfections. However, in actual practice the test data ex- 
hibit scatter with reliable values of lower magnitude. The 
eccentricities which precipitate buckling of unpressurized 
cylinders again are present. 

It is interesting to see what happens to the buckling load 
as we increase the value of the internal pressure. At first 
there is a steady sharp rise in the buckling stress. Then when 
some value of pressure is reached, the buckling stress no longer 
increases but remains approximately constant. Thereafter 
an actual decrease in allowable buckling stress takes place as 
we increase pressure. This latter behavior may be explained 
as follows. It is Known that eccentricity is detrimental in 
the compressive buckling of cylinders. When pressure is 
high, the structure actually has a built-in eccentricity due to 
the bulging of the cylinder. This action is prevented at the 
supports. Eventually, as one approaches the yield strain of 
the material in the hoop direction the compressive strength of 
the cylinder should drop rapidly. The structure then be- 
haves like an anisotropic shell. For these reasons it is be- 
lieved that one should use the strain parameter (pR/Et)/en, in 
plotting the design curves. Here ¢,, is the tensile yield strain 
of the material in the hoop direction. 


Fig. 4b Compressive buckling stress as a function of pressure 
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Fig. 5 Bending buckling stress as a function of pressure hoop 
strain, L/R = 1 


Some investigators have used the increase in the buckling 
stress coefficient over the unpressurized one as the ordinate in 
plotting the design curves. It is not considered of advantage 
to plot this change in buckling coefficient as ordinate. A more 
realistic and direct picture is obtained if one uses the actual 
buckling coefficient. In this way a family of curves is ob- 
tained as shown in Figs. 4a and 4b. However, one must 
distinguish between values of L/R. The weakening effect of 
length carries over for this loading condition, since the effect 
of eecentricity is one of the “hidden” parameters. However, 
percentagewise the effect of length is less than for the unpres- 
surized loading. For example, from Figs. 4a and 4b, for 
R/t of 2500, the decrease in allowable stress in going from 
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Fig. 6a Torsion buckling stress as a function of pressure hoop 
strain for cylinders with simply supported edges, L/R = 1 
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L/R of 1 to L/R of 6 is about 17 per cent. The corresponding 
decrease for the same ¢ylinders when unpressurized is about 
37 per cent according to Fig. 1b. 

For conventional materials, the tensile yield strain (o,,/E) + 
0.002 is about 0.007. Therefore this value has been used 
for oy, in preparing the curves in Fig. 4. None of the avail- 
able test values of pR/Et exceeded 0.0025 [(pR/Et)/0.007 ~ 
0.35]. More tests should be carried out for values of pR/Et 
greater than 0.0025 to establish the correct shapes of the 
curves in Figs, 4a and 4b. It is of interest to note that pre- 
vious investigators suggest the use of the parameter (p/E) X 
(R/t)? which stems from the theoretical considerations. The 
following example is given to illustrate the shortcomings of 
the use of that parameter. Suppose we have two cylinders of 
the same material, say 2024 ST aluminum, one with R/t = 
500, the other with R/t = 5000. If we let p = 100 psi for 
the first cylinder and p = 1 psi for the second, (p/E) (R/t)? 
is 2.5 for both cylinders. However, the first cylinder pre- 
sumably has failed, since pR/t = 50,000 psi. Meantime the 
second cylinder is working at a hoop stress of only 5000 psi. 

The curves for R/t = 5000 in Figs. 4a and 4b have been 
drawn to fit the general trend of the family of curves in that 
figure. Test data are needed to verify this trend. 


> 


Combined Bending and Internal Pressure 


This loading condition occurs when a side load is applied to a 
missile tank. The test data are limited (25). Again there 
is a large scatter. The basie behavior of the structure under 
pressure is similar to that of the axially compressed pressur- 
ized cylinder. There is a rapid rise in allowable stress with a 
relatively small amount of internal pressure, followed by a 
leveling off. However, the maximum allowable stress for this 
loading condition exceeds the one for the axially compressed 
pressurized cylinder. The physical reasons for the increased 
buckling strength are as follows. First, as already explained 
for the unpressurized case, the maximum stress is localized. 
Secondly, a cylinder in bending “flattens”? or changes shape 
in cross section. This type of deformation is opposite in di- 
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Fig. 6b Torsion buckling stress as a function of pressure hoop 


strain for cylinders with simply supported edges, L/R = 6 
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rection to the bulging action due to the internal pressure, at 
least in the regions of maximum stress. Therefore, eccen- 
tricity due to internal pressure is lessened by this counter- 
action. The rise in allowable stress coefficient recommended 
for design is shown in Fig. 5 for R/t = 3000. The rise for 
R/t = 1000 is not clearly defined because of lack of test data. 
It appears to be at a smaller rate than for R/t = 3000. 
Additional test data are needed to confirm this trend. 


Combined Torsion and Internal Pressure 


There are enough test data to draw design curves as shown 


in Figs. 6. The increase due to internal pressure for this 
loading condition is very high and does not appear to approach 
a limit with increasing pressure as in the previous cases of 
compression and internal pressure and bending and internal 
pressure. However, the allowable stresses appear to ap- 
proach a constant envelope curve independent of R/t for 
large internal pressure. Here again the parameter used is 
pR/Et. More test data are needed for pR/Et > 0.002 to 
establish the true nature of these curves. The present curves 
should not be extrapolated. 


Combined Transverse Shear and Internal Pressure 


No test data are available for this loading condition to the 
knowledge of the author. In light of the behavior of unpres- 
surized cylinders under transverse shear, it might be recom- 
mended that the allowable stress be taken tentatively as 1.25 
times the buckling stress for the loading of combined torsion 
and internal pressure discussed previously. This factor may 
be all right for cylinders with low internal pressure and low 
R/t. However, in actuality, the factor could turn out to be 
unity for cylinders with high internal pressure and high R/t. 


Multiple Combined Loadings 


Test data here are meager for cases involving internal pres- 
sure. Some data exist for combined compression, torsion 
and internal pressure, and combined compression, bending 
and internal pressure. The scatter of the test results is about 
the same as the scatter one obtains for unpressurized cylin- 
ders under combined loadings. Furthermore, the shapes of 
the interaction curves are similar to the curves representing 
the interaction of loading conditions for unpressurized cylin- 
ders. The curve shown in Fig. 7, which is reeommended for 
design, is therefore valid for cylinders under combined loads 
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Fig. 7 Buckling for combined compression, bending, shear, 
torsion and internal pressure 


with and without internal pressure. Notice that the curve 
has been extended beyond both axes. The upper left-hand 
section represents combined tension and bending with little 
torsion. The lower right-hand section represents combined 
tension and torsion with little bending. 


Combined Compression and External Pressure 
According to theory (4), a linear interaction equation is 

suggested for this case. 

Combined Torsion and External Pressure 


Use the equation (51): R, + Rr! = 1. 


Areas for Further Test _ 


The design data presented in the previous sections have 
been based on existing known test data plus extrapolation in 
areas where no data were available. Therefore further tests 
must be performed for the following conditions in order to 
insure the reliability of design curves: 

1 Unpressurized cylinders in axial compression with 
R/t > 1500 and L/R > 3. 

2 Unpressurized cylinders in bending with R/t > 1500 
and L/R > 3. 

3 Pressurized cylinders in axial compression with pR/Et > 
0.0025; also R/t > 3000. 

4 Pressurized cylinders in bending with pR/Et > 0.0025; 
also R/t > 3000. 

5 Pressurized cylinders in torsion with pR/Et > 0.002. 

6 Pressurized cylinders in transverse shear. 

7 Pressurized cylinders under combined loadings such as 
compression, bending and shear; tension, bending and shear; 
bending and torsion. 

Any additional test information on the buckling of cylin- 
drical shells also would be worth while. For example, all pre- 
vious data are based on test specimens made up of sheets 
longitudinally spliced. Tests on cylinders made by circum- 
ferential splicing of sheets would be of interest. 


E = Young’s modulus 

L = length of cylinder a @ 

M = applied bending moment , a 

P = applied axial load 

p = pressure 

R = radius of cylinder + Li 

S = applied transverse shear force 

T = applied Torque ‘5 

= thickness of cylinder wall 

Z = length parameter, 0.95 (L/R)? R/t Oo 

e = strain 7 

o = axial stress 

+ = shear stress 

Subscripts op 

b = bending 

c = compression 

cr = critical 

hy = hoop yield 

r = reduced 

t = tangent 
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HE REALIZATION and employment of space flights 

are based on a theoretical analysis and on numerical cal- 
culations of the equations of motion for vehicles in cosmic 
space. The basic requirements for the powering and guid- 
ance of the launching rocket, as well as the optimal and the 
acceptable conditions for the takeoff, are determined by com- 
putation. 

As the initial coordinate system for describing the travel of 
cosmic vehicles, we may take a Cartesian coordinate system 
with the origin at the Earth’s center, which is in translational 
motion relative to the stars. In many practical problems 
spherical coordinate systems have to be used, rigidly fastened 
with the Earth, with the origin at the center of the Earth and 
at various points on the Earth’s surface. The movement of 
bodies relative to the moon or other planets is also to be con- 
sidered. 

Into the differential equations of motion of celestial ballis- 
tics describing free flight in cosmic space, we must introduce 
only the integration forces which are determined by Newton’s 
universal gravitation law. In solving problems of hitting the 
moon and flight around the moon, we must consider the body’s 
movement in a certain gravitational field produced by the 
sun, the moon and the Earth, taking into account the oblate- 
ness of the Earth. 

In groping for the optimal conditions of a rocket’s takeoff 
for the moon, approximation methods may be used, permit- 
ting us to consider the body’s Keplerian motion relative to 
Earth when the distance from the moon is more than 66,000 
km, and as Keplerian motion relative to the moon when this 
distance is less than 66,000 km. 

For given points in the northern hemisphere, the optimal 
conditions for a takeoff in a first approximation when only 
the Earth’s gravitational forces are considered may be ob- 
tained in the following way. 

First, let us consider the collision case. Let A be the posi- 
tion of the point of cutoff of the rocket, B the position of the 
moon’s center at the moment of collision, and O the Earth’s 
center (see Fig. 1). Any orbit corresponding to a flight from 
point A to point B lies in a plane determined by the three 
points AOB. For a full determination of the orbit and the 
magnitude of the initial velocity V, it is sufficient to specify 
the inclination angle & to the horizon of the initial velocity V. 

From the solution of the two-body problem, we easily find 
the functional relation 


3 


where V, is the parabolic velocity. The angle between the 
directions OA and OB, denoted as ®, is called the angular dis- 
tance. Point A, which is rigidly tied to the Earth, wasn 
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a latitude tne of the Earth as sit rotates. The angle b changes 
within certain limits, which are seen in Fig. 1. 

Since the-ratio r/R is small and practically constant, ac- 
cording to function [1] the ratio V?/V,? depends essentially on 
Jand . In Fig. 2, this dependence on #, with = constant, 
is graphically represented by the continuous lines; the dotted 
lines represent the typical form of the rocket characteristics 
for various constant values of the rocket weight at the mo- 
ment of cutoff. 

It can be seen from Fig. 2 that the optimal starting time, at 
which the useful weight is the greatest, corresponds to the 
maximum possible value of the angular distance ®, this condi- 
tion determining the most advantageous starting time for 
the given point A. Evidently, in this case the plane of the 
orbit is normal to the plane of the Equator. The considered 
maximum value of the angle @ differs for various positions of 
the moon in its orbit around the Earth. Evidently, for points 
in the northern hemisphere, the maximum value of ® cor- 
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responds to the lowest position of the moon below the equa- 
torial plane. The above conditions define the optimal time 
within a moon month, as well as the launching moment 
within a 24-hr period. 

The previous conclusions were made without taking into 
account the rotation of the Earth and the advantage of using 
an additional velocity component of the rocket due to the 
rotation of the Earth. If the rotation of the Earth is taken 
into account, the inclination of the optimal orbit to the equa- 
torial plane is reduced. Some corrections must also be made 
on account of the influence of the moon, the sun, the Earth’s 
oblateness and the characteristics of the powered part of the 
flight. 

An additional substantial increase in the payload can be 
secured by the choice of the launching point on the Earth’s 
surface. These most favorable points lie near the Equator. 

For the Soviet cosmic rockets, a version with an inclination 
of the initial part of the orbit to the equatorial plane amount- 
ing to 65 deg was decided, which is near to the most advan- 
tageous one. 

The basic calculations of the equations of motion were 
made on high speed computers. Optimal trajectories were 
determined, and a set of trajectories closest to the optimum 
region was studied. 

The choice of any one orbit may be begun by stating the 
total specific energy of the missile at the moment of cutoff. 
The value of the required total energy is determined by the 
mission of the flight and by a need for insuring favorable ob- 
servation conditions for watching the missile from the terri- 
tory of the USSR at the time of its hitting or passing near the 
moon. This also determines the duration of flight until the 
moon is reached. An analysis of the optimal conditions 
makes it evident that the time of flight from the Earth to the 
moon area must be approximately 0.5, 1.5, 2.5, 3.5 or 4.5 
days. 

After the required total energy has been determined, it is 
essential to have information about the possible payload 
weight and the necessary launching conditions for several 
days close to the most favorable day. The latter corre- 
sponds, for the northern latitudes, to the least inclination of 
the moon at the moment of the closest approach. This 
problem is answered through ballistic computations in which 
the rocket’s characteristics are taken into consideration. 

The computations show that a miss of 2 to 3 days on both 
sides of the optimal launching time keeps the useful weight 
loss within acceptable limits. Therefore only several days 
within each moon month appear to be the proper time for 
launching. A considerable deviation of the launching day 
from the optimal one causes major reduction of the possible 
payload weight. 

For a detailed study of the influence of initial parameters 
on the trajectory characteristics, it appeared helpful to es- 
tablish some essential properties of a set of trajectories, of 
which the central trajectory hits the center of the moon. 
The following properties of a set of trajectories, determined as 
in the foregoing, are sufficiently accurate; the use of these 
properties greatly facilitates the solution of many practical 
problems. 

1 Near the moon, in a limited area around the central 
trajectory which may be considered as straight line, the 
nearby trajectories in the motion relative to the moon over 
distances up to 20,000 km, form a set of trajectories that have 
a rotational symmetry about the central axial path. 

2 In the vicinity of the moon, the trajectories are approxi- 
mately conic sections which lie in meridional planes. 

3 The various orbits of this family may be obtained for 
small deviations from the principal values of the six independ- 
ent parameters. The coordinates of the rocket and the 
components of its absolute velocity at the moment of cutoff 
may be considered as the six parameters. 

4 A change in the increment of any individual parameter 
or a proportional change in any combination of the six param- 
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any 


eters leads to various orbits of the set, all of which lie in the 
same meridional plane. 
5 A general nondimensional variable \ can be introduced 
for all meridional planes, equal to the ratio 


where 


Ag = increment of any parameter & at the end of launching 
Aé* = increment corresponding to the given minimal distance 
of the perturbed orbit from the center of the moon 


Consequently, a universal relation 
Pmin = f(A) 


is valid in any meridione! plane, for all various orbits. Here 
Pmin is the minimal distance of the path from the center of 
the moon. 

Using the properties previously indicated, the six-parameter 
set is converted to a two-parameter one, which may be mapped 
uniquely on an arbitrary plane 7, in which the points are de- 
termined by the following polar coordinates: \X—the radius 
and w—the polar angle. w stands for the angle defining the 
meridional plane of the orbits. 

This circumstance allows us to examine and describe the 
aggregate of perturbed movements through a certain clearly 
observable procedure. 

It is convenient to draw in the plane 7 isolines of the vari- 
ous numerical characteristics of the trajectories. This is 
particularly important for examining a flight-around orbit, and 
for finding initial data corresponding to the orbits that have 
the required and the most advantageous characteristics. 
An example of two families of isolines is presented in Fig. 3; 
these were used to determine the orbit of the flight—a round 
trip of the Automatic Interplanetary Station (AIS). 

The families of isolines of maximum and minimum dis- 
tances from the object to the Earth during the first revolution 
after approaching the moon have been plotted on this figure. 
The initial energy for the chosen set of trajectories corre- 
sponds, in the absence of the moon’s influence, to an apogee of 
550,000 km from the Earth’s center. For those trajectories 


of the set that lie to the left of the isoline of 550,000-km apo- 
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Cz = projection on the axis of vernal equinox xX 


Cz = projection on the Earth’s rotation axis. 


Complete 
geocentric For the geocentric angular momentum 
energy Cy C= + Cy? + 
Prior to passing the —0.68 km?2/sec? 55 X 10%km?2/see —30 x 10° 21 X 10° 67 X 103 
influence sphere of _) 
the moon me 
After passing the in- —0.74 -177X 108 65x10? 32X10% X 108 


Cy = projection on the axis perpendicular to the x axis in the equatorial plane. 


gee, deceleration occurs and the apogee decreases; to the right 
of this line, acceleration takes place, the apogee increases, and 
the energy of the object with respect to the Earth rises after it 
approaches the moon. The corresponding quantitative effects 
can well be seen on the diagram. 

The values of pmin = 7900 and and w, which were real- 
ized during the motion of the AIS, have been determined by 
computation. By plotting these data (Fig. 3), we obtain an 
apogee of 480,000 km, and a perigee height on the first revolu- 
tion equal to 47,500 km. 

The perturbing influence of the moon on the motion of the 
automatic station is very great; it can be characterized quan- 
titatively by the data found in Table 1. 

It is to be noted that the magnitude of the vector of angular 
momentum increases approximately threefold and its direc- 
tion is almost reversed. 

The z plane simplifies the analysis of the scattering effect of 
the initial parameters on the orbit characteristics. The in- 
troduction of the 7 plane allows us to combine the problems 
of hitting the moon and of flights close to the moon. The 
problem of determining flight-around trajectories with pre- 
scribed characteristics is reduced to that of hitting a given 
point in the 7 plane, corresponding to the required values of 
the orbit characteristics. 

Flight around the moon can be realized along trajectories 
of different types. 

The trajectory of a remote flight around the moon (at 
distances of about 40,000—-100,000 km from the moon, with 
only slight influence of the latter) is near to an ellipse with a 
focus at the center of the Earth. Having this trajectory and 
starting from the northern hemisphere, the return path to- 
ward the Earth after passing the moon will be below the 
Equator. This makes it impossible to observe the object 
from the territory of the USSR when it is close to the Earth. 
Besides, for these trajectories the perigee distance from the 
center of the Earth is less than the radius of the Earth, and 
therefore the object enters the terrestrial atmosphere on the 
first revolution. Hence, its lifetime is short. 


It was found to be of advantage to realize a trajectory of 
another type, with the object passing near the moon at a dis- 
tance of about 5,000-20,000 km. In this case, the moon 
causes a big perturbation, as a result of which one can ob- 
tain a trajectory returning to the neighborhood of the Earth 
and subsequently have a satellite with a high perigee of about 
40,000 km and an apogee of about 500,000 km, the return to 
Earth taking place over high latitudes of the northern hemi- 
sphere, which creates very favorable conditions for observa- 
tion and reception of radio information. In addition, such 
trajectories are very convenient for solving the problem of ob- 
ject orientation and photographing at distances ranging be- 
tween 40,000 and 150,000 km. Calculations have shown 
that within this range of distances, which is covered in about 
12 hr, the direction of the object-to-moon vector in absolute 
space remains practically unchanged. This is an advantage 
for normal operation of the orientation system and for pho- 
tography. This feature of the orbit allows switching on of the 
orientation system and the subsequent photographing opera- 
tions by means of a program device operating from a preset 
timing mechanism. 

In order to photograph the invisible side of the moon, it is 
necessary to launch the rocket on days near to the new moon. 
In that case, the invisible side of the moon is illuminated by 
the sun. The combination of favorable energy conditions 
with those of the illumination of the moon determines the 
months and the days which are most suitable for obtaining 
photographs of the back side of the moon. 

The foregoing considerations served as a basis for choosing 
the flight-around trajectory for the third cosmic rocket, mak- 
ing use of the strong influence of the moon for the best solu- 
tion of the problem in question. 

The actual data concerning the orbits of the three Soviet 
cosmic rockets are given in the following and summarized in 
Table 2. 

1 The last stage of the first cosmic rocket weighed 1472 kg 
without the fuel; the weight of the container, scientific ap- 
paratus and energy sources amounted to 361.3 kg. The 
rocket was launched toward the moon on Jan. 2, 1959; it 


Table Soviet cosmic rockets 
Name of Launching Weight of last stage Weight of scientific apparatus ; ' 
No. rocket date without fuel, kg and power supply, kg re Purpose a 
1 Lunik I Jan. 2, 1959 1472 361.3 artificial planet of solar system 
2 Lunik II Sept. 12, 1959 1511 390.2 to moon, Sept. 14, 1959 
Pe wre 0 hr 02 min 24 sec 
Moscow time 
3 Lunik III Oct.4,1959 ay, 1553 435 (weight of AIS—278.5) AIS photography 
Tes of invisible side of moon 
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passed north of the moon at a distance of about 5000 km from 
its surface. Flight time from the Earth to the proximity of 
the moon was 1.5 days. 

It was in the launching of this rocket that parabolicspeed was 
first exceeded, and consequently, after passing near the moon, 
the rocket continued to move away from the Earth, leaving the 
sphere of the Earth’s gravitational field and becoming an arti- 
ficial satellite of the sun. 

Its orbit around the sun has the following characteristics: 


Inclination of the orbit plane to the ecliptic plane i = 1 deg. 
Minimum distance from the sun 146 X 106 km. 

Maximum distance from the sun 197 X 106 km. 
Corresponding period 7 = 450 days. 

- Minimum distance from the sun was reached in the middle 
of January 1959. The shortest distance from the orbit 
of the artificial planet to the orbit of Mars is 15 X 106 km. 


Radio communication with the first cosmic rocket was main- 
tained up to distances of 400,000-500,000 km. 

2 The second cosmic rocket was launched with the aim of 
hitting the moon. The weight of the last stage without fuel 
was 1511 kg; the weight of the container with the scientific 
apparatus, 390.2 kg. The motion diagram is presented in Fig. 


Plane of lunar orbit 
Pattern of flight to the moon 


Fig. 4 


The following data characterize the flight orbit of this 
rocket: The trajectory lies practically in one plane making an 
angle of 65 deg with the Equator. Flight from the Earth to 
the moon lasted about 1.5 days. The initial velocity in free 
flight exceeded the local parabolic speed. 

With respect to the precision of initial data, a flight time of 
about 1.5 days proved to be more advantageous than schemes 
of hitting after longer flight times (2.5, 3.5 or 4.5 days), for 
which lesser initial energy is required. For a flight to the 
moon within about 0.5 day, too great energy consumption is 
required. 


/ Comet of encounter 
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The minimum geocentric speed on the orbit slightly ex- 
ceeded 2 km per sec. The velocity of impact with the 
moon’s surface was 3.3 km per sec. The velocity vector 
made an angle of about 60 deg with the moon’s surface. 

The point of impact on the moon is situated about 800 km to 
the north of the center of the visible lunar disk, south of the 
craters of Archimedes, Aristilles and Autolycus. The point 
of impact on the lunar map is shown on Fig. 5. The cor- 
responding pmin for the prolonged trajectory is 500 km. 

The moment of impact was registered from an abrupt 
ceasing of radio signals on Sept. 14, 1959, at 0 hr 2 min 24 
sec a.m. Moscow time. 

This moment of impact agrees with the data of trajectory 
measurements which were carried out both during the flight 
and during the time interval directly preceding the hitting 
moment. 

A projection of the trajectory of the container of the second 
cosmic rocket on the Earth’s surface is given in Fig. 6. 

3 The third cosmic rocket was launched on Oct. 4, 1959, 
for a flight around the moon. The weight of the last stage 
without the fuel was 1553 kg. The weight of the Automatic 
Interplanetary Station (AIS) with the apparatus and the 
energy sources, installed in the last stage of the missile, 
amounted to 435 kg. 

During the flight of the AIS from the Earth to the moon the 
inclination of the orbit to the equatorial plane was equal to 
55 deg. After the perturbation by the moon, further move- 
ment under the influence of the Earth’s gravitation took place 
along a near elliptic orbit with an inclination to the Equator 
of about 80 deg. The calculation of further movement shows 
that the sun and the moon influence the AIS orbit so that the 
inclination of the orbit changes irregularly and decreases. 
On the tenth revolution, the inclination amounts to 48 deg. 
On the eleventh revolution, the inclination of the orbit under 
the influence of the moon increases again to 57 deg. It is 
noteworthy that the minimal distance from the Earth de- 
creases from one revolution to another due to the influence of 
the sun and the moon. A calculation shows that after com- 
pleting the eleventh revolution at the end of March, the AIS 
will enter the Earth’s atmosphere in the northern hemisphere, 
thus ceasing to exist. 

This circumstance is due to the orbit form and the char- 
acter of its location relative to the Earth and the sun. This 
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effect, though unexpected at the first glance, depends only 
upon Newtonian forces. It is evident that such effects should 
be taken into account during theoretical analysis of problems 
concerning the structure of planet systems and the properties 
of the orbits of different planets and their satellites in the 
solar system. As a result of perturbation due to the sun, an 
evolution takes place which may lead to the collision of a 
satellite with its principal planet; therefore only satellites 
with orbits of certain types can “survive” for a long time. 

Fig. 7 gives projections of the first revolution of the AIS 
orbit on the Equator plane and on the plane which is perpen- 
dicular to the direction from the Earth’s center to the vernal 
equinoctial point. 

Fig. 8 gives the projection of one revolution of the AIS on 
the Earth’s surface. The projection line shows loops when 
passing near the Earth’s rotation axis both in approaching 
the Earth and in moving away from it. 

Fig. 9 shows the orbit projection on the moon’s surface 
(continuous line); points are marked signifying the minimal 
distance to the moon’s center and the photographing time. 
A grid of selenographic coordinates is plotted on the moon’s 
surface. The limits of the area which has been photographed 
by the AIS are marked with dotted line; the hatched part 
lies in shadow, the white part is sunlit. 

The conditions for photographing were favorable with re- 
spect to the illumination of the invisible side of the moon, 
and an area comprising two thirds of the surface of the invis- 
ible side was photographed. An unfavorable factor in the 
process was direct illumination of the surface and, as a result, 
the absence of shadows from the accidents of the moon’s sur- 
face. The picture was obtained only as a result of varying 
reflectivity of its different parts. 

Before photographing, the orientation system was auto- 
matically switched on, turning one end of the axis of the AIS 
to the sun and the other end with the objectives of the cameras 
toward the moon. After stabilizing, the hatch covering the 
illuminator and protecting the transmitters of the moon 
orientation and the cameras was opened. After fixing the moon 
direction with a special photocell device, the light seeking de- 
vices reacting to the sunlight were automatically switched 
off; subsequently, the accuracy of the orientation of camera 
axis to the center of the moon was increased. After complet- 
ing the orientation to the center of the moon, the process of 
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photographing began at 06.30 hr, Moscow time, on Oct. 7, 
the AIS being then at a distance of 65,500 km from the moon’s 
center. To insure favorable temperature conditions inside 
the AIS after the photographing process, the station was forced 
to rotate around its transverse axis with an angular velocity 
of 2 deg per sec. Thereupon the film was automatically de- 
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veloped, fixed and dried. The picture of the moon was trans- 
mitted to the Earth by means of a special television system. 
Figs. 10, 11 and 12 give samples of photographs processed 
to date. 
The pictures transmitted to the Earth from aboard the 
AIS show not only the side of the moon invisible from 
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Fig. 8 Projection of the orbit of the AIS on the Earth. Numbers on the diagram correspond to the successive positions of the pro- 
jection of the orbit of the AIS on the Earth’s surface 


1-16 hr, 8/4/59, 131 thous km from Earth; 2- 

4-16 hr, 8/5/59, 266 thous km from Earth; 5- 

7-16 hr, 8/6/59, 361 thous km from Earth; 8- 
10-16 hr, 8/7/59, 410 thous km from Earth; 1l- 
13-16 hr, 8/8/59, 443 thous km from Earth; 14- 
16-16 hr, 8/9/59, 464 thous km from Earth; 17- 
19-16 hr, 8/10/59, 473 thous km from Earth;  20- 
22-16 hr, 8/11/59, 472 thous km from Earth; 23- 
25-16 hr, 8/12/59, 459 thous km from Earth; 26-— 
28-16 hr, 8/13/59, 436 thous km from Earth; 29- 
31-16 hr, 8/14/59, 400 thous km from Earth;  32- 
34-16 hr, 8/15/59, 349 thous km from Earth;  35- 
37-16 hr, 8/16/59, 281 thous km from Earth;  38- 
40-16 hr, 8/17/59, 186 thous km from Earth;  41- 


43-16 hr, 8/18/59, 


52 thous km from Earth; 44- 


46— 0 hr, 8/20/59, 189 thous km from Earth 


Hemisphere not seen from Earth Hemisphere facing Zartn 
Gy \ \ 


Ohr, 8/5/59, 183 thous km from Earth; 3- 6hr, 8/5/59, 216 thous km from Earth; 
Ohr, 8/6/59, 300 thous km from Earth; 6- 6hr, 8/6/59, 324 thous km from Earth; 
Ohr, 8/7/59, 279 thous km from Earth; 9- 6hr, 8/7/59, 392 thous km from Earth; 
Ohr, 8/8/59, 423 thous km from Earth; 12- 6hr, 8/8/59, 432 thous km from Earth; 
Ohr, 8/9/59, 452 thous km from Earth; 15- 6hr, 8/9/59, 457 thous km from Earth; 
0 hr, 8/10/59, 468 thous km from Earth; 18- 6 hr, 8/10/59, 471 thous km from Earth; 
0 hr, 8/11/59, 474 thous km from Earth; 21- 6 hr, 8/11/59, 474 thous km from Earth; 
0 hr, 8/12/59, 469 thous km from Earth; 24- 6 hr, 8/12/59, 466 thous km from Earth; 
0 hr, 8/13/59, 453 thous km from Earth; 27- 6 hr, 8/13/59, 447 thous km from Earth; 
0 hr, 8/14/59, 425 thous km from Earth; 30- 6 hr, 8/14/59, 416 thous km from Earth; 
0 hr, 8/15/59, 384 thous km from Earth; 33- 6 hr, 8/15/59, 372 thous km from Earth; 
0 hr, 8/16/59, 329 thous km from Earth;  36- 6 hr, 8/16/59, 312 thous km from Earth; 
0 hr, 8/17/59, 253 thous km from Earth; 39- 6 hr, 8/17/59, 230 thous km from Earth; 
0 hr, 8/18/59, 148 thous km from Earth; 42- 6 hr, 8/18/59, 114 thous km from Earth; 
0 hr, 8/19/59, 55 thous km from Earth; 45-16 hr, 8/19/59, 150 thous km from Earth; 
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the Earth, but a small portion of the moon visible from the 
Earth as well. This permitted connecting heretofore un- 
known objects on the moon’s surface with those which are 
familiar, and to determine their coordinates on the moon. 
Precise shapes have now been obtained for the objects situ- 
ated at the edge of the visible disk, little known so far because 
of a distorted perspective. In other cases the objects were 
known only partially, and the pictures received allowed us to 
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get more precise shapes for the parts situated on the visible 
side and to discover the shapes of their unknown continua- 
tions on the invisible side. 

A mountain range, named the Sovietsky Mountain Range, 
and the presence of seas and craters, were discovered on the 
moon’s invisible side. One of the seas has been named the 
Sea of Mechta, in memory of the first cosmic rocket launched 
on Jan. 2, 1959. 

Figs. 10 and 11 show the original photographs, from which 
the obvious effects of radio noise have been cleared. The 
picture shown in Fig. 12 was obtained as a result of prelimi- 
nary treatment of a number of photographs, representing a 
first attempt at a composite based on details revealed in sev- 
eral photographs. For a complete description of the invis- 
ible side of the moon, the photographing must be continued. 
It can be fairly expected that further pictures, with side illu- 
mination, will reveal more unknown formations on the other 
side of the moon. 

Weare of the opinion that this experiment of photographing 
in space and subsequently transmitting the pictures to Earth 
is interesting not only from the point of view of studying the 
nature of the back side of the moon, but is also of essential 
importance as a way of testing new methods in modern experi- 
mental astronomy. 


Fig. 12 Physical features of the lunar hemisphere invisible 
from the Earth, determined after preliminary processing of 
the photographs obtained from the AIS: 1 A large crater 
sea 300 km in diameter—the Sea of Moscow. 2 The Gulf 
of Astronauts of the Sea of Moscow. 3 Continuation of 
Mare Australe on the back side of the moon. 4 A crater with 
a central peak—Tsiolkovsky. 5 A crater with a central peak— 
Lomonosov. 6 Jolio-Curie Crater. 7 Sovietsky Mountain 
Range. 8 The Sea of Dream. 

Physical features of the visible side of the moon: I Mare 
Humboldt; II Mare Crisium; III Mare Marginis; IV Mare 
Underum; V Mare Smythii; VI Mare Foecunditatis; VII 
Mare Australe. 

The continuous line across the chart is the lunar equator. The 
broken line indicates the border between the hemispheres visible and 
invisible from the Earth. Features circled with a continuous line 
have been definitely established as a result of preliminary examina- 
tion; features circled with a broken line are to be further studied 
with regard to their shape; features circled with a dotted line are 
being studied with regard to their classification; with regard to the 
other parts, the photographs are being studied further. 
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Some Problems of Providing for 
Scientific Research on Rockets 
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ODERN rocketry provides the scientist with a powerful 
tool for probing mysteries of nature heretofore inacces- 
sible to man. Starting with vertical rocket launching, fol- 
lowed by placing artificial satellites in orbit with the aid of 
rockets, and finally launching cosmic rockets for the investiga- 
tion of cosmic space in the vicinity of the moon. Soviet 
scientists obtained within slightly more than a decade sub- 
stantially new data on the properties of the Earth’s upper 
atmosphere, on the phenomena that take place there, on the 
composition of cosmic radiation, and on the properties of inter- 
planetary matter. As is known, the research performed to 
date culminated in the first photographs of the dark side of 
the moon. 


Soviet scientists owe this success both to the vigorous de- 
velopment of rocketry proper, and to rapid progress in many 
fields of science (radio engineering, electronics, physics and 
chemistry of metals, etc.). Planned and systematic work was 
carried out in the Soviet Union to perfect the rockets them- 
selves as well as the scientific apparatus and procedures for 
the scientific research, along with broadening the research 
programs. In addition, each successive experiment not only 
enriched science, but provided additional experience for fur- 
ther improvements to the necessary technical apparatus, and 
raised new technical problems in turn. At the same time, 
each experiment showed the way to new techniques and con- 
structions, and helped eliminate numerous difficulties. The 
paper by Prof. Krasovskii discloses many scientific results 
obtained in the course of this research. My own task is to 
show several particular solutions that made possible research 
with the aid of rockets. 


Even the very first rocket investigations of the Earth’s 
atmosphere demonstrated an urgent need of isolating the 
measuring instruments from the rocket carrier, so as to get 
rid of disturbances and inaccuracies caused by rocket motion. 
It was also necessary to eliminate mutual influences of some 
instruments on other, i.e., to separate the devices from one 
another in space. Therefore, in almost all experiments, the 
scientific apparatus was installed in compartments isolated 
from the rocket itself. 


Apart from isolating the payload compartment with the 
scientific apparatus, certain containers were ejected by 
special mortars, symmetrically located along the longitudinal 
axis of the rocket. Fig. 1 shows the placement of the mortars, 
which were equipped with pneumatic propulsion devices 
actuated in flight at an instant of time predetermined by the 
program. Such devices were used on certain geophysical 
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rockets during the IGY. The containers acquired an addi- 
tional speed relative to the rocket, at a small angle with its 
longitudinal] axis, and were thus kept at a satisfactory distance 
from the rocket and from each other. 


Fig. 2 shows the flight trajectory of another type of instru- 
ment container, which was to be parachuted safely to the 
Earth (for example, carrying devices for obtaining air samples 
at certain altitudes). A general view of such a container is 


shown in Fig. 3. It consists of three parts (shown in Fig. 4). 
The head part is pressurized and contains electric batteries, 
devices for registering the readings of the manometers, a 
camera for photographing these readings, illumination for the 
scales of the gages, the program mechanism, amplifiers for 
the ion manometers, and cameras for taking pictures of clouds 
and the surface of the Earth and neighboring space. 


Screens 


| Payload container 
Mortar 
3 Pneumatic cylinder 


ARS JouRNAL 


: 
| 
; 
> 
~ 
2 
‘ 
CORES, 
> 
ae 
| 
| 
= 
— 


200 
180 


160 
| H-2 xn 
120 — 


are provided on the outside for registering meteor particles, 
There is also a bumper with a corrugated metal device for 
softening the landing impact. A fanning system of rods also 
serves the same purpose, and sometimes keeps the landing 
container in a vertical position. The middle compartment of 
the container is a latticed shell, containing glass bulbs for 
air sampling. The third part contains the parachute system. 


Fig. 3 1—Parachute compartment, 2—central compartment, 3—instrument compartment, 4 
6—storage battery compartment, 


5—nose compartment, 
ites, 


partment, 
8—manometers, 
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Upon landing, the deformation of the head bumper generates 
a signal which uncouples the parachute, so as to prevent the 
wind from dragging the container along the ground. 

This type of container is naturally not the only one. 
Mortars were also used to eject containers that did not need 
to be landed, for example, those carrying the telemetering 
system for the instrument readings. The principal advantage 
of the demonstrated system is therefore the possibility of 
performing different research programs with standard rockets 
without rebuilding the latter. All that need change is the 
equipment in the containers. 


A second example of the problems encountered in the re- 
search concerns the orientation of the objects ejected in space, 
a task which also necessitated a number of investigations. 

For example, one type of geophysical rocket carried a 
spherical container with a motor driven gyroscope wheel in 
its equatorial plane (Fig. 5). This gyroscope was accelerated 
directly before the launching of the rocket. After the con- 
tainer was separated from the rocket, the gyroscope axis 
remained almost vertical. Two photoelectric transducers 
were located in the equatorial plane 180 deg apart to provide a 
rough orientation. Depending on their position relative to 
the sun’s rays, a signal was generated either to decelerate or 
to accelerate the gyroscope wheel, thereby tilting the container 
about the vertical axis in a direction necessary to bring the 
sun’s ray to the null meridian between these transducers, 
where the instrument for precision orientation was installed. 
The latter disconnects the rough-orientation transducers at 
that instant and takes over control of the acceleration or 
deceleration of the wheel, thus aiming the null meridian 
directly at the sun (Fig. 6). In all the experiments performed, 
good orientation of the container was achieved until the latter 
entered the dense layers of the atmosphere. 

Another method of orienting an object was used in the last 
launching of the lunar rocket to photograph the moon’s 
surface from an automatic interplanetary station (AIS). 
(See Fig. 7.) The jet principle was used to rotate the station: 
Compiessed air escaped from miniature nozzles, two pairs of 
which served to produce rotation about the longitudinal axis, 
and one about the transverse axes. Signals to start the gas 
flow from the various nozzles were generated by a special 
logical computer, which developed the optimum nozzle 
operating program from information on the direction of the 
sun’s rays, obtained from the photoelectric transducers. 

The flight trajectory of the AIS shows the successive steps 
in the operation of the orientation system (Fig. 8): I—Start 
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of orientation; IIl—steadying the station (elimination of 
rotation); I[1[—station oriented toward the sun; IV—moon- 
orientation transducers are switched on automatically a cer- 
tain time interval after the system is steadied; these trans- 
ducers are located on the opposite side behind the window, 
the angle between the longitudinal axis aimed toward the sun 
and the orientation of the moon being 13 deg. Precise orien- 
tation toward the moon occurs in the section IV-V. Pho- 
tography begins in position V. After the film advance is ter- 
minated, the window shutter is closed. In position VI, a 
pulse is produced to rotate the station to ensure further nor- 
mal operation of the solar battery. In position V, one can 
see the field of view of the lenses of the two cameras. 

As a third example of a solution of an individual technical 
problem raised by this research, we describe the heat control 
system used in the AIS. 

The heating conditions of an AIS moving in orbit are 
determined by the joint action of the external and internal 
heat sources. The only external heat source is the sun. 
The internal heat is generated by the equipment. For normal 
operation of the AIS equipment, its temperature must be 
kept within fixed limits. 

The AIS is a cylindrical body with two spherical end covers. 
Located inside the AIS is scientific, television and radio 
apparatus. The body of the station carries solar batteries, 
and on its cylindrical part is mounted a radiation surface 
(Fig. 9, 1) which is covered by rotating blinds (Fig. 9, 6) in 
the form of perforated cylindrical sleeves. The heat from 
the instruments is transferred to the shell through the frame 
by radiation and by forced circulation of the gas. The latter 
is produced by a fan (Fig. 9, 4) through an air duct (Fig. 9, 2). 

The blinds open whenever the gas temperature rises above 
+25 C and close when the gas temperature drops below a 
permissible value. The use of blinds permits measurement 
of the optical coefficients of the outer surface of the AIS and 
thus regulates its heat conditions. 

The nonradiating portion of the shell surface (Fig. 9, 3) is 
covered with screens, which serve to equalize the temperature 
field over the shell. 


The power supply to the research and radio apparatus was 
always a troublesome matter and the subject of special care. 
The use of solar batteries on satellites and cosmic rockets has 
substantially eased the situation. Let us illustrate the general 
solution of the problems of power supply to the AIS. 

The AIS power system contains autonomous chemical 
power sources to feed apparatus with a short operating time, 
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as well as a centralized buffer block of chemical batteries 
(Fig. 10, 3). The spent energy of the buffer cells is replen- 
ished by the solar batteries (Fig. 10, 1), which are charged 
between broadcasts. During the time of the broadcast, the 
apparatus is fed from both the buffer and the solar batteries, 
via a distributing network. 
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Of particular interest in this system is the choice of shape 
of solar battery, which in turn influenced the overall view of 


the AIS. The shape selected insures high power stability in 


unoriented flight. By placing the elements of the solar bat- 
tery on the surfaces of a 12-face prism and two truncated 
pyramids, the summary effective area of the entire battery, 


Electric motor 


= which is proportional to the output current, changes very 
little with changing direction of illumination. 

The figure shows the variation of power with the angle 
of incidence of the rays. In the formula of Fig. 11, S; is the 
illuminated area, a; the angle of incidence of the ray on the 
face, and n; the number of faces. 


Finally, we dwell briefly on the organization of the physio- 
logical research, aimed at investigating the influence of various 
factors on living organisms during rocket flight. As is known, 
many flights of animals (principally dogs) were organized in 
the Soviet Union. These flights, too, necessitated the solu- 
tion of many technical problems. 

Let us describe a pressurized cabin for dogs and its equip- 
ment, as used in one type of geophysical rocket, designed for 
altitudes up to 475 km. The cabin permitted maintenance 
of normal pressure, gaseous air composition and relative 
humidity. It contained two removable trays for holding 
the animals, a regeneration plant and recording apparatus. 
The cabin was placed in the parachute equipped compartment 
of the rocket. The cabin was shaped nearly like a truncated 
cone, with bases in the form of spherical segments. The 
volume of the cabin was 460 liters. The inner wall of the 


a cartridge and control manometers (Fig. 12). In addition, 
to keep the air pressure in the cabin constant, a pressure 
regulator was used. In subsequent experiments, high activity 
chemical compounds were used to regenerate the air in the 
cabin. These were located in a specially constructed regener- 
ation cartridge (Fig. 13). Draft was produced in the car- 
tridge by an electric fan, which delivered 20 liters per min. 
Two cartridges were installed in the pressurized cabin, making 
the total volume of air 40 liters per min. With such a regen- 
erating plant, the concentration of carbon dioxide in the 
cabin did not exceed 0.7 per cent, and the oxygen concentra- 
tion remained practically constant. 

To register the physiological functions, a set of medical 
apparatus was used, consisting of an automatic pressure 
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recorder, an amplifier block, a matching network and a set of 
electrodes and transducers. 

Records were made of the arterial pressure of the oscilla- 
tions of the carotid (by means of a piezoelectric transducer), 
and of the depth and frequency of respiration. The respira- 
tion was recorded with a transducer comprising a rubber 
plate covered with a thin strain gage, the resistance of which 
increased with stretching of the rubber. The amplitudes of 
the records of breathing motion were regulated with a potenti- 
ometer. 

The electrocardiograms of the animals were taken with the 
aid of a block consisting of four amplifiers feeding the vibra- 
tors of the automatic recorders. 

In addition to recording the physiological functions, motion 
pictures were taken of the behavior of the dogs. A special 
electrically driven camera (fed with 26-30 v d-c) operating at 
24 frames per sec was used. The camera cassettes held 300m 
of film, enough for 11 or 12 min. The camera was provided 
with a short-focus lens of 20-28 mm focal length. The 
camera was mounted in the upper part of the cabin with the 
lens pointing downward. To increase the optical distance 
between the lens and the dogs, the pictures were taken 
through mirrors (Fig. 14). 

The animals were illuminated with two 60-w bulbs. A 
watch with a second hand was included in the field of view 
of the camera. 

The equipment installed in the cabin for a prolonged flight 
of an animal in a satellite—that of the historic Laika—was 
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HE ACHIEVEMENTS of the Soviet Union in rocketry 

made it possible to put into practice some new and impor- 
tant scientific explorations of the upper atmosphere and cosmic 
space. Soviet engineers and scientists created gigantic 
Sputniks and cosmic rockets. Geophysicists and astrophysi- 
cists placed their instruments on them to research dif- 
ferent properties of the upper atmosphere and cosmic space, 
hard electromagnetic and corpuscular radiation of the sun, 
interplanetary dust und gaseous medium, magnetic fields, 
cosmic rays, planets and their satellites, and different proc- 
esses around them. Before the launching of Sputniks, Soviet 
scientists were experienced in the exploration of the upper 
atmosphere with the help of geophysical rockets. 

The launching of the first Soviet Sputnik marked a new 
era—an intensive assault aimed at complete mastery of the 
secrets of the upper atmosphere and interplanetary space and 
at the creation of continuously working laboratories beyond 


‘This paper was presented, along with the papers by Sedov 
and Blagonravov, at the ARS Annual Meeting, Nov. 17, 1959, 
Washington, D. C., in response to an official invitation of the 
Society. 

2 Chief, Department for Research in Upper Atmospheric 
Physics, Institute of Atmospheric Physics. 
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Results of Scientific Investigations Made f = 
by Soviet Sputniks and Cosmic Rockets 


more complicated. Particularly stringent requirements were 
imposed here on the pressurization of the cabin, which was 
in the form of a cylinder 650 mm in diameter and 800 mm 
long. 

In addition to the equipment for air regeneration and 
temperature control, it became necessary to provide for 
automatic feeding, sanitation equipment, and telemetry of 
the air pressure, of the temperature at various points of the 
cabin, of electrocardiogram records, breathing depth, oscilla- 
tions of the cartoid, air pressure in the sphygmomanometer 
bag, and finally, the movements of the animal, which had 
limited freedom to lie down, stand and sit. 

The data received made it possible to establish the variation 
of the physiological functions of the animal both during the 
active ascent of the rocket and during the orbiting of the 
satellite. 

A discussion of these results is beyond the scope of this 
paper, particularly since a translation of an article on this 
subject was published in the October 1959 issue of ARS 
JourNAL. However, I mention this example to illustrate 
one of the numerous sets of problems that must find technical 
solution in the organization of scientific research. Often 
the success or failure of an investigation depends on the 
solution of similar problems, which at first glance appear to 
be of secondary importance. 

The success of the research now carried out in the USSR 
is evidence of an immense effort exerted toward solving 
numerous varied technical problems. 
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the Earth’s atmosphere. Every new Sputnik was more per- 
fect than its predecessor and enriched our knowledge about the 
upper atmosphere and cosmic space. 

The use of cosmic rockets is a further improvement in this 
respect. With their help, the scientists managed to get a 
fuller representation of huge processes in space surrounding 
the Earth, for the first time reaching another cosmic body— 
the moon—and getting information not known before. 

The purpose of the present paper is to describe the results of 
scientific researches carried out by the Soviet Sputniks and 
cosmic rockets. Large groups of scientists, engineers and 
workers of many scientific institutions, design offices and in- 
dustrial enterprises were engaged in this great research. All 
this became possible due to the high level of development 
of Soviet economics and industry. Valuable information 
about the upper atmosphere and cosmic space has been ac- 
cumulated up to the present moment. 

I have been instructed by all my colleagues to make a 
survey of the results of the subjects already treated. Much of 
it is still being treated and carefully studied. 

Let us start with a few words about scientific research of 
the upper atmosphere. Before the launching of the first 
Soviet Sputnik, the upper atmosphere beginning with the 
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200-km level was thought to be isothermal and to have a 
temperature of less than 1000 K, and, consequently, the scale 
height was thought small. However, by deceleration of the 
Sputniks, it was found that at a height of 220 km, the density 
of the atmosphere was equal to 2.5 X 10- gm cm=*,_ It 
exceeded what had been expected, and this result essentially 
changed the conceptions of density and temperature of the 


Fig. 1 Photo of an artificial sodium cloud at a distance of 

152,000 km from Earth during the flight of the second cosmic 

rocket. The pictures were taken by an electron telescope with 

an exposure “5 sec. The increase in diameter of the cloud 

corresponds to the sequence in time. The interval between the 
photos is 10 to 20 sec 
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Fig. 2 Survey of different values of density of the upper at- 
mosphere. The left curve represents data obtained by means 
of manometers; the right, by other means 


upper atmosphere accepted before. By the change in the 
period of revolution of Sputniks around the Earth and by the 
change in the position of the perigee, it was found that, at the 
220-km level, the density of the atmosphere was greater by 
day than at night and greater in the polar than in the equa- 
torial regions. 

To determine the exact coordinates of the cosmic rockets, a 
sodium cloud was tried. By means of a sodium vapor cloud 
the coordinates of the first and second cosmic rockets at a dis- 
tance of 100,000—-150,000 km were determined with an ac- 
curacy of a few angular seconds. The time of cloud formation 
was registered with an accuracy of some seconds in time. By 
way of preparation for this experiment, a sodium cloud was 
brought to being at the level of 430 km in twilight. Fig. 1 
shows different stages in the formation of sodium clouds. 
By data on diffusion of this cloud, it was found that 1.8 X 
10° particles cm~ are present at this height. The compari- 
son of values obtained with those from the retardation of the 
Soviet Sputniks showed that the estimated scale height at the 
430-km level is 70 km. This points to the fact that the tem- 
perature of the upper atmosphere within the experimental 
region is 1600 K. For sucha great scale height the density of 
air at a given level is much higher than was supposed previ- 
ously. 

Ionization and magnetic manometers were also mounted on 
the third Sputnik and their use allowed data about the density 
of the upper atmosphere to be obtained which agree with a 
new representation of the properties of the upper atmosphere 
within the error. 

Fig. 2 shows the indicated properties of the upper atmos- 
phere. 

The third Soviet Sputnik carried a spectrometer for deter- 
mining the nature of the atmospheric ions with a mass number 
of 6 up to 50 units. Before mounting this instrument on the 
Sputnik, it was tested on high altitude rockets. The ion of 
nitric oxide was found to be the most frequent up to the height 
of 250 km. The existence of this ion up to this height does 
not mean that it forms as a result of ionization of the neutral 
molecule of nitric oxide. As was shown by a number of 
theoretical and laboratory investigations, the ionized mole- 
cules of the nitric oxide appear in the upper layers of the at- 
mosphere as a result of the reaction of an oxygen ion with a 


~ neutral molecule of nitrogen or in the reaction of the oxygen 
~ atom with an ionized molecule of nitrogen. The nitrogen ions 


found show that nitrogen molecules are present at altitudes up 
to 250 km. However, the exploration with the help of the 
third Sputnik revealed that there are only atomic ions of 
oxygen and nitrogen higher than 250 km; the ions of atomic 
oxygen are dominant, whereas the ions of atomic nitrogen con- 
stitute only several per cent of the oxygen ions. This is evi- 


_ dence that the upper atmosphere above 250 km has an atomic 


structure. 

The third Sputnik and the cosmic rockets carried special 
ion traps which make it possible to estimate the ion density in 
the upper atmosphere and cosmic space. If a certain assump- 
tion is made as to the mechanism of charging the rocket and 
characteristics of the ion current depending on various dif- 
ferences of potentials between the electrodes of the ion trap, 
then the material obtained shows that the bodies are being 


_ charged in the ionosphere up to a negative potential of several 


volts. It has been shown that the estimated electron density 
is several million electrons cm~* up to the height of 475 km, 


_ and that this is in good agreement with earlier investigations 


by ultra-short wave dispersion interferometers on highaltitude 
rockets. It has also been found that there is a slow decrease 
in the density of ionized particles with height, and that above 
the maximum of the F,; region, there are twice or three times 
as many electrons as below. 

On the basis of the data about the electron density from the 
study of reception of radio signals from Sputniks, similar con- 
clusions have been drawn. All the results are compared in 
Fig. 3. 
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All the cosmie rockets carried ion traps. They helped to 
determine the currents between the electrodes of ion traps at 
different distances from the Earth’s center in cosmic space and 
near the moon. The preliminary treatment of the results 
allows us to assume that the density of ionized particles in the 
outer atmosphere at a distance of some radii from the Earth 
is of the order of several thousand ion pairs per cm? or less. 
Somewhat higher currents between the electrodes of an ion 
trap were observed near the moon. This requires further 
study. 

Sputniks and cosmic rockets were used for the study of 
cosmic rays, using Geiger counters and scintillation counters 
combined with different absorbers. Fig. 4 illustrates the re- 
sults of these investigations in comparison with the altitude 
dependence of the intensity of cosmic rays available since 
1949. It concerns the particles capable of penetrating 10 gm 
em~*. The instruments of the second Sputnik made it pos- 
sible to determine the intensity of cosmic rays as a function of 
the geomagnetic latitude. The results of these investigations 
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Fig. 4 The intensity of hard cosmic rays, penetrating 10 gm- 
cm~*, as a function of height. The solid curve represents the 
results obtained with rockets, the dashed curve above 400-km 
data obtained by the second Sputnik. The scale along the ordi- 
nate axis is in relative units 
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are given in Fig. 5. On Nov. 7, 1957, at 04:40 (Moscow 
time), an abnormally large increase in the intensity of cosmic 
rays was registered at high geomagnetic latitudes due to an 
external radiation zone (see Fig. 6). It was the first record of 
hard particles within an outer belt. There were no other 
explorations of this kind before at such a height. 

Still more valuable information was obtained by the third 
Sputnik on regularity of change of intensity of pulse counting 
due to electrons of about 100 kev. Fig. 7 is a chart of data 
for May 15-22, 1958. A rapid growth of intensity of count- 
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Fig. 5 Latitude relationship of hard cosmic rays, obtained with 
the help of the second Sputnik 
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discovered by the second Sputnik, Nov. 7, 1957. The dashed 
curve reflects the normal variation of the intensity, and the solid 
curve observed values 
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Fig. 7 Latitude distribution of intensity of 100-kev electrons. 
Black circles represent the position of the Sputnik when the 
instruments went off scale on the ascending loops and crosses the 
position of the Sputnik when they went back on scale in the 
descending loops. The dashed line is a geomagnetic parallel 
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ing is denoted by circles. Crosses represent the emergence of 
the Sputnik from the polar zone, where a higher intensity was 
observed. The crosses southward point to the relationship 
of the effect with height, since the crosses correspond to the 
backward loops, i.e., the greater height of flight. It was found 
that the zone of high intensity is asymmetrical with respect to 
the Earth’s magnetic dipole. For illustration, the geomag- 
netic parallel is denoted with a dashed line. The component of 
100 kev does not have a constant intensity. The analysis of 
the signals from the third Sputnik in the vicinity of South 
America showed the existence of a region of higher intensity of 
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Fig. 8 Intensity of ionization in a scintillation counter as a 
function of a geomagnetic latitude and height above Earth’s sur- 
face. The ordinate axis denotes energy in ev/sec deposited 
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Fig. 9 Intensity of hard radiation as a function of distance from 
Earth during the flight of the first cosmic rocket. N—number of 
pulses registered by the counter; E-—energy deposited in a 
crystal of sodium iodide. Roman numerals denote data: 
I, II and III—rate of counting with a threshold of 45 kev, 450 
kev and 4.5 mev, respectively, with rate of counting referring to 
the area unit of the cross section of a crystal (19 cm?); IV—total 
ionization (total power delivered by a crystal per sec); V—read- 
ings of the Geiger counter. Number of impulses refers to the 
area unit of the cross section of counters (4 and 15 cm?). For 
great intensity, the data of a small counter are used. The posi- 
tion of the orbit of the first cosmic rocket with respect to radiation 
belts is shown schematically at upper right. Time corresponding 
to given points of the orbit is shown along the orbit. Distri- 
bution of the intensity of corpuscles at different distances from 
the Earth is shown schematically above the trajectory 


ionizing particles above the Equator. The recorded particles 
in this region were harder than in the vicinity of polar regions. 

Fig. 8 is a record obtained in the southern hemisphere in th« 
range of heights from 1600 down to 1100 km. It is seen from 
the figure that, as the Sputnik moves toward the Equator, 
the intensity strongly increases although the height of th« 
flight essentially decreases. The numerous analyses show 
that this equatorial zone of the higher intensity of cosmic rays 
extends up to 45 deg of geomagnetic latitude. Hard protons 
are, obviously, of importance here. In the material published 
it was indicated that, no matter what the mechanism of form- 
ing the particles in the equatorial zone is, it is apparent that 
the factor of accumulation plays the main role in this effect. 
The concentration of particles in the equatorial zone, where 
they are able to oscillate at a great height above the Earth for 
a long time, presents convincing evidence. 

The investigation of hard ionized particles showed the es- 
sential qualitative and quantitative difference between hard 
radiation near the equatorial zone and that of the polar zone. 
All these results were reported at the 5th meeting of CSAGI in 
Moscow in August 1958, and at the International Conference 
on Peaceful Use of Atomic Energy at Geneva in September 
1958. 

Further investigation of hard particles revealed the distri- 
bution of intensity of particles of different energy at various 
distances from the Earth in the vicinity of the moon. Fig. 9 
shows the results obtained by the first cosmic rocket launched 
on Jan. 2, 1959. They enabled estimation of the type of 
energetic spectrum of electrons in the range of 20 to 100 kev. 
If the integral energetic spectrum of electrons is represented 
in the form 


N(> 


then y ~ 5. The spectrum is notably harder (y ~ 8) for the 
sections of smaller intensity at the outer edge of the zone. 
This points to an increase in the number of particles as the 
energy decreases. The maximum of intensity of the outer 
zone of the increased radiation was found to be 1} radii 
nearer to the Earth at the moment of exploration by the second 
cosmic rocket (Sept. 12, 1959) than by the first. By applying 
absorbers of different materials (lead and copper), the hard 
electron component was observed to contain in this zone a 
significant quantity of particles of energy greater than 0.5 
mev. 

It was thought before the launching of Sputniks that illu- 
mination, heating and ionization of the upper atmosphere are 
mainly due to the hard electromagnetic radiation of the sun. 
It was assumed only for the polar regions that the solar cor- 
puscules (charged particles: Protons, a particles and elec- 
trons) may penetrate into the atmosphere while geomagnetic 
disturbances and aurorae are occurring. The study of the 
spectra of these aurorae showed that they are frequently 
caused by protons with a sizable range of speeds. However, 
no hydrogen emission was observed in many cases. Here 
again the aurorae were presumably explained by not very 
hard electrons of up to hundreds and thousands of electron 
volts penetrating into the atmosphere. 

An attempt was made to discover these not very hard elec- 
trons by the third Sputnik. Two very thin fluorescent 
screens with aluminum foils of different thicknesses were used 
as indicators. The emission of the fluorescent screens was 
received by a photoelement. The photocurrent was then 
amplified. An electric signal was transmitted to a memory 
system and then, by radiotelemetry, transmitted to Earth. 
By using metal foils of different thickness it was possible to 
estimate both the intensity and the energy of the electrons 
that produced the greatest fluorescence on the screen. The 
peculiarity of the instrument was in its reaction practically 
only to electrons, without recording protons and photons of 
the same energies. As a result, powerful fluxes of electrons 
of about 10 kev were recorded at heights up to 1900 km over 
the South Pacific. The intensity of these fluxes was very high, 
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Fig. 10 Reading of the indicator of electrons with an aluminum 
foil 0.8 < 10~* gm cm~? as a function of an angle between the 
axis of the indicator and a magnetic field line. The scale of the 
indicator Ko.s is proportional to the logarithm of the current of 

electrons le * 


h 
and in the majority of cases the instruments appeared to go off 
scale, since such high intensity was not expected. The effec- 
tive energy of electrons changed during an experiment. The 
intensity of fluxes showed a tendency to decrease as the effec- 
tive energy of the particles increased. The effective energy 
of electrons was observed to be less in the polar regions. The 
energy flux of the electrons studied exceeded 100 erg cm~? 
sec! at the height up to 1900 km at the moment of going off 
scale. If the electron fluxes of the indicated intensity had 
penetrated into the lower layers of the atmosphere, i.e., the 
F region of the ionosphere, they could not have remained un- 
noticed, since they would essentially have increased ionization 
of the upper atmosphere and led to the appearance of aurorae. 
Since no such phenomena were observed, the recorded fluxes 
were explained to be electrons oscillating along the magnetic 
field lines (July 1958, Fifth Meeting of CSAGI). 

Thus, information was obtained about the accumulation of 
electons of about 10 kev at great heights up to 1900 km. All 
other explorations known to us refer to much lower altitudes 
and cannot be an indication of the existence of a magnetic 
trap around the Earth. 
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Fig. 11 Effective energy of electrons as a function of the angle 
between the indicator’s axis and a magnetic field line, with data 
in this diagram corresponding to Fig. 10. The data represent 
the position of the Sputnik near the apogee within the range of 
height of 1720-1880 km. This section is divided into three 
intervals. The closed and open circles denote the consequent 
data, referring to these intervals. The symbol ! denotes off- 
"} scale values which actually can be a little greater 


The majority of not very hard electrons move along the 
directions normal to the magnetic field lines. The electron 
flux is greater downward than upward. This can be seen in 
Fig. 10. In Fig. 11, the dependence of effective energy of 
the discovered electrons on the direction of the electron motion 
with respect to the magnetic field lines, is illustrated. Since 
the angle subtended by the indicators was } steradian, the 
diagram shows that the opposite motion of electrons from the 
Earth near the magnetic field lines is practically absent. 
The electrons moving to the Earth were observed even under 
small angles to the magnetic field lines. This indicates that 
the particles penetrating to the lower layers of the atmosphere 
appear as a result of some processes at heights exceeding 1900 
km. It has been found that the energy of electron fluxes able 
to reach the F layer of the ionosphere without reflection can 
reach values of about 1 erg cm~? sec™!. Partial increases of 
intensity of corpuscles were registered above the Pacific 
even at —4 deg of geomagnetic latitude at the height of 
~1500 km, as was already reported in the material published. 

It is no exaggeration to say that the powerful flux of the 
electrons discovered in the upper atmosphere is of extreme 
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importance in understanding many important processes, al- 
though it creates no threat for astronautics. It is of interest 
to note that sizable intensities of such electrons begin to ap- 
pear at that geomagnetic latitude where, in the F region of the 
ionosphere, there had also been observed previously the in- 
crease of ionization but which could not be explained by the 
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Fig. 12 Readings of the Cherenkov counter as a function of the 

distance from Earth’s surface. Data denoted by a dashed line 

refer to the flight of the first cosmic rocket, and by the solid line 

to the second cosmic rocket. Scale of the ordinate axis is in 
angular units 
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Fig. 13 Readings of nuclei with a charge equal to 2 or more asa 
function of distance from Earth. Scale along the ordinate axis is 
in relative units 
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Fig. 14 Changes in intensity of the magnetic field during the 
third Sputnik’s flight above the East Siberian world anomaly. 
Time in sec is measured along the abscissa axis; height and 
angle of the indicators’ position along the ordinate axis. Num- 
bers show the data referring to: 

1 Field directly measured by the magnetometer above USSR 
territory. 

2 Field measured by a magnetometer after elimination of 
the influence of the magnetic deviation caused by magnetic 
parts of the instrument. 

3 Field at Earth’s surface along the trajectory of the flight, 
taken from the magnetic map. 

4 and 5 Readings of the orientation transmitters, showing a 
precession with a period of 136 sec and 
axis with a rate of 0.36 deg/sec. yO Ss 

6 Change of flight height. 
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effect of hard electromagnetic radiation of the sun. Some in- 
homogeneities in the ionization of the upper atmosphere and 
the system of ring currents around the Earth can be explained 
by the existence of such corpuscles. The dependence of the 
temperature variations and the density of the upper atmos- 
phere upon solar activity becomes more understandable now, 
since the corpuscles, the intensity of which is governed by 
the solar activity, are the sources of heating. The heat flux 
brought by these electrons can well account for the increase of 
temperature with height and latitude, which could not be 
ascribed to the hard electromagnetic radiation of the sun 
alone. 

The problem of great accumulation of hard particles in the 
upper atmosphere is now being widely discussed among geo- 
physicists and astrophysicists, and there are different view- 
points on the problem. 

Harder corpuscles discovered in the upper atmosphere are 
of no geophysical importance. They cannot essentially change 
the state of the upper atmosphere since their energy content is 
small. However, these hard corpuscles, as well as x-radia- 
tion arising in the Earth’s atmosphere and in the body of 
rockets and Sputniks while they are irradiated by these elec- 
trons of energies of tens of kilo electron volts, are very un- 
pleasant for astronautics, since the dose of dangerous radia- 
tion may reach tens of roentgens per hour. It is impossible 
to neglect this radiation, which requires that special measures 
be taken to defend cosmic travelers or photoemulsions in the 
apparatus. The investigations of hard radiation at great 
heights show that flight from the Earth of future astronauts 
will be less dangerous if the trajectory of the flight passes 
through the polar regions having the least concentration of 
hard particles. The best time for the flight from the Earth 
is, obviously, the time immediately after strong geomagnetic 
disturbances, after which the content of hard particles in the 
upper atmosphere greatly decreases. 

The cosmic rockets carried Cherenkov counters of charged 
particles of a very great energy. The application of discrim- 
inators permitted registration of nuclei with charges greater 
than or equal to 2, 5 and 15. Apart from this, thanks 
to a photoeffect in the photomultiplier, x-rays of energy of 
tens and hundreds of kilo electron volts were registered. 
Since this radiation was mainly due to bombardment of the 
rocket bodies by hard electrons, then the registration of the 
sum photocurrent reflected the intensity of these particles. 
Figs. 12 and 13 show the results of these investigations. The 
intensity variation of x-radiation reflects the passage of indi- 
cators through the radiation belt around the Earth. In con- 
formity with the experiments described in the foregoing, the 
maximum of the registered intensity during the flight of the 
second cosmic rocket was approximately 13 radii lower than 
in the flight of the first cosmic rocket. The fluxes of positive 
nuclei with charges greater than or equal to 2, 5 and 15, 
changed at great distances from the Earth, are in a ratio of 
250:20:1, respectively. It is of interest that counting of 
nuclei with a charge equal to 2 or larger is of a greater value in 
the center of a radiation belt around the Earth. This cir- 
cumstance is now being carefully studied. The existence of 
a particles in the radiative zone is supposed to be one of the 
possibilities. 

The third Sputnik carried a magnetometer to investigate 
the Earth’s magnetic field at great heights. This instrument 
determined the orientation of the Sputnik in the space around 
the Earth. Anomalies within the geomagnetic field were dis- 
covered by it. Many of them are of short duration and, ob- 
viously, are related to the ionospheric currents which cause 
magnetic disturbances. It is worth noting also that the char- 
acter of height changes of intensity of abnormal and normal 
fields above the East Siberian world anomaly is the same, and 
this may be an indication of the fact that its sources are very 
deep. (See Fig. 14.) It allows one essentially to specify the 
known theories of a geomagnetic field and bring them into 
accord with the actual data. 
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It is of special interest to investigate the magnetic field of 
the Earth with the help of a magnetometer carried by the first 
and second cosmic rockets. Fig. 15 illustrates the expected 
change of the intensity of the magnetic field and its actual 
value during the flights of these rockets. So, for the first 
time, it was possible to register the magnetic field up to a dis- 
tance of 60,000 km from the Earth, created by charged parti- 
cles in the magnetic trap of the Earth. The result is in good 
agreement with the existence of not very hard electrons dis- 
covered around the Earth and described in the foregoing. 

Attention should be drawn to the fact that the actual values 
of the magnetic field, given in Fig. 15, fluctuate. The devia- 
tions are much greater than those expected from the error in 
instruments. Thus, there are grounds to suspect that changes 
of intensity of the magnetic field near the Earth, estimated 
to be up to some tens and even a hundred of gammas, really 
exist. These fluctuations are being carefully studied. There 
are grounds to suppose that they are caused by short-period 
variations of the geomagnetic field, which were already 
thought by us to be exciters of the discovered corpuscles. 
Qur attention is attracted by the closeness of the values of 
density of such a variable magnetic field and fluxes of elec- 
trons with an energy of about 10 kev at the moments of their 
maximum intensity. It is quite possible that these forms of 
magnetic and corpuscular energy are in balance or close to it. 

To register micrometeors, the third Sputnik and the cos- 
mic rockets carried piezoelectric elements. These elements 
registered the value proportional to impulse caused by a mi- 
crometeor’s shock, which is, in turn, proportional to its en- 
ergy. Table 1 gives information about the number of regis- 
tered micrometeors and about the flux of matter brought by 
them into the atmosphere. 

Exploration of the moon is most remarkable. It was es- 
tablished by the second cosmic rocket that the magnetic field 
close to the surface of the moon does not exceed 50-100 y. 
The information about the magnetic field in the vicinity of 
the moon allows us to state that the effective magnetization 
of the moon is not higher than 5 per cent compared to the ef- 
fective magnetization of the Earth. This is new and very 
important information for the conception of magnetism of 
cosmic bodies. We observed no inverse cube variations of 
the magnetic field near the moon. It is of interest to note, 

however, that fluctuations of the magnetic field similar to 
those discovered at the distance of some radii from the Earth 
have also been registered. Further investigations will reveal 
whether these fluctuations really exist and what characteris- 
tics they have. Moreover, if they reflect real values of the 
magnetic field, frozen into the interplanetary gas, then it is of 
great interest from the point of view of astrophysics, since it 
is the first direct indication of the interplanetary magnetic 
field. It has also been discovered that, within the error of 


Note: Material is treated partially. 


Table 1 Intensity of micrometeors 


—_ mass of registere¢ ntensity of a flux matter in tons 
Means of - > micrometeors in gm, of micrometeors for the whole 
exploration Date = 4 X 10% cm m~* sec“! globe per 24 hr 
Third Sputnik May 15,1957 = 2x 107% 5 to 10 5 to 10 x 10° 
May 16-17, 1957 (8.10-* to 30.10-°) 5 xX 10-3 5 xX 108 
May 19-26, 1957 <10-* <10? 
First cosmic rocket 10~° <2 <10 
Second cosmic rocket 10-* x 10-* ~102 
107 to 20. 10-%) 
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Fig. 15 Intensity of the magnetic field 7 as a function of the 

distance from the center of Earth. Curve 1 represents data ob- 

tained during the flight of the first cosmic rocket, curve 2 of the 

second cosmic rocket. Curve 3 shows the variation of the ex- 

pected variation of magnetic field with distance, taking into 
account the dipole and quadrupole terms 


measurements made close to the moon, there is no increase of 
intensity of hard corpuscular radiation as compared with that 
at the distance of some lunar radii. Cherenkov counters of 
hard charged particles also did not discover near the moon 
any deviations from that registered at a distance of several 
radii. The first photos of the unseen side of the moon have 
been obtained, and they will contribute to knowledge of the 
nature of this cosmic body. 

As a result of the explorations described, our knowledge of 
the upper atmosphere and cosmic space have been greatly ex- 
tended. New phenomena have been discovered and quanti- 
tative data obtained about those assumed before. A fascinat- 
ing perspective of the exploration of the outer atmosphere 
and cosmic space at ever increasing heights, and of planets 
and their satellites by new and perfect means, now faces us. 
All this is a basis for new theoretical investigations. 

A lively discussion is now taking place on the problems of 
physics of the upper atmosphere and cosmic space, and the 
not-too-distant future will give the answers to all these ques- 
tions. 
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tors, in addition to structures. 
requirements is emphasized. 


HE BOLD venture of manned flight into space presents 

sizable challenges in the conception and design of space 
vehicle structures. For the space vehicle must be for its 
occupants a world away from their natural world. It must 
provide a life-supporting atmosphere and protect them against 
meteoroid bombardment and the hazards of radiation; it 
must maintain comfortable cabin temperature levels in the 
presence of prolonged exposure to solar radiation; it must be 
a storehouse of food and water to nourish men, and of fuel 
to power machines; it must provide for the elimination of 
their waste products and for the expulsion of the large amounts 
of heat given off by electrical, electronic and mechanical com- 
ponents; it must provide an environment conducive to the 
maintenance of healthy mental and physical dispositions 
under otherwise unstimulating conditions. This means that 
for long duration flights there must be ample room for mov- 
ing about, adequate noise suppression and reasonably com- 
fortable furnishings. In addition, the space vehicle may well 
have to supply its own pseudogravitational field for the 
comfort and perhaps the safety of occupants during prolonged 
flights. Moreover in certain specific aspects, it must be more 
than our natural world, since it must survive the inferno of 
entry into planetary atmospheres while preventing the 
occupants from being baked to crisps! Finally, it must 
perform all of these functions having the minimum possible 
weight. 

It is evident, therefore, that the consideration of space 
vehicle structures demands the compatible integration of 
many and varied technical disciplines. From the earliest 
stages of conception of the space vehicle, structural considera- 
tions, such as materials and construction operating under a 
wide range of temperatures, must be treated intimately 
with aerodynamic, thermodynamic, guidance, control and 
human factors considerations. For instance, a winged entry 
glide vehicle with sharp leading edges, although desirable 
from the viewpoint of aerodynamic efficiency, might be in- 
tolerable because of the enormous rate of aerodynamic heating 
over a sharp edge. Even though such a design were feasible 
from the thermodynamic consideration of available cooling 
schemes, it might not be feasible structurally because of 
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Structural Considerations of 
Manned Space Vehicles 


This paper considers the structural requirements for a manned space vehicle. 
leading to an optimum vehicle considering the requirements of thermodynamics and human fac- 

The necessity of early and thorough integration of these differen: 
Several space vehicle configurations are presented including a ballis- 
tic and a glide type, in order to exemplify some of the various structural problems encountered in 
each type during the phases of spaceflight, entry into an atmosphere and landing. 
aerodynamic heating and loading, meteoroid impact and other problems. 
tions suitable for optimum space structures are discussed sufficiently to demonstrate available 
choices and indicate areas of required development. 
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excessive thermal stresses and strains. Increasing the lead- 
ing edge radius, however, though structurally feasible, 
might result in such intolerable aerodynamic losses, that : 
complete redesign of the basic shape might be required. 
Such an integrated approach is, of course, not a novel concept 
by any means. We are all familiar with the design cari- 
catures which result when they are made from an isolated view- 
point. But the emphasis here is that this integration of 
disciplines must be considerably more intimate than in 
ordinary aircraft practice and timed sufficiently early to 
achieve optimum performance with minimum weight. 

The problem of structural design of space vehicles, there- 
fore, is more involved than merely formulating adequate 
structure on the basis of specified and uncoupled loading 
conditions; it is rather one part of an evolutionary process 
whereby the optimum vehicle design is attained through 
suitable tradeoffs with many and often widely different 
technical viewpoints. 

In this paper some of the structural considerations of 
manned space flight such as atmospheric entry, heat protec- 
tion schemes, meteoroid protection and landing will be dis- 
cussed. Present state of the art and development needs will 
be indicated where pertinent. 


Entry Into the Atmosphere 


We first examine the entry phase of space vehicles into a 
planetary atmosphere, say the Earth’s, since it is in this 
phase that the various disciplines of aerodynamics, thermo- 
dynamics, control and structures must be intimately bound 
together in order to achieve optimum design. The vehicle 
structural configurations meeting the requirements of the 
entry phase will then be considered in relation to the way 
they meet some of the other requirements stated previously. 

Since we are concerned with manned vehicles we shall 
confine our attention to blunt nosed ballistic vehicles having 
low weight to drag ratios (W/CpA) and glide vehicles having 
moderate lift to drag ratio (L/D), all of whose initial entry 
path angles are low. These conditions tend to minimize the 
severity of deceleration loading and aerodynamic heating. 

In general a ballistic vehicle will experience higher levels of 
deceleration than a glide vehicle. Typical deceleration peaks 
attainable in manned ballistic and glide vehicles (L/D~1) 
are 10g and 1g respectively without the aid of drag devices. 
These decelerations do not present difficult structural prob- 
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lems and are within the human tolerance limitations. A 
comparison of the heating rate histories of the various types 
of vehicles is given in Fig. 1. The heating rate curve cor- 
responding to the pure ballistic type is for the stagnation 
point, since this represents the major heating problem for a 
blunt nosed vehicle. The curve given for the glide type 
(made available to the author by Flathers) corresponds to the — 
side walls, since due to the very large proportion of side wall 
area to stagnation areas, this is representative of the major 
heating problem affecting the structure. The representative 
maximum heat flux of the ballistic type is considerably larger 
than that of the glide type, whereas the total heating time 
of the former is much smaller than that of the latter. The 


significance of this comparison is that blunt nosed ballistic 


vehicles must be designed to withstand relatively high rates _ 
of heat input with less total heat transferred to the cabin. 

Beryllium heat sink or ablation appear to be low weight 

solutions for handling the heat flux shown in Fig. 1. The glide — 
types, on the other hand, must be designed to handle a rela- 
tively large total heat input, necessitating perhaps radiation- 
type cooling for the greater portions of the area. The 
stagnation regions like leading edges of wings, however, re- 
quire protection schemes capable of handling much higher | 
leat fluxes. For instance a 3-in. radius leading edge on the - 
glide vehicle considered in Fig. 1 would have a heating rate 


approximately 10 times the peak side wall value for a sub- _ 


stantially long time. 
To facilitate discussion of some of the structural aspects of — 
« ballistic-type entry vehicle, let us consider the hypothetical | 
examples shown in Figs. 2 and 4. In both cases the vehicles 
are assumed to have a large radius blunt frontal surface 
called the heat shield. In Fig. 2 the shield considered is a 
heat sink type of beryllium. Beryllium, because of its rela-_ 
tively large heat capacity per unit weight and high strength, 
is outstanding as a heat sink despite its low ductility at room | 
temperature and around 1200 F (1)? and the current problems 
of its toxicity and manufacturing of large one piece units. 
Since a heat sink sustains a high overall temperature and 
therefore undergoes a great deal of expansion, the manner 
in which it is mounted to the remaining structure is of con- 
siderable importance. A possible configuration (Fig. 2a) 
has the shield fastened to the supporting structure by means 
of a number of bolts which are free to move radially with the 
expanding shield. A layer of resilient material and insulation 
provides a distributed contact between the shield and sup- 
porting structure. Designing in this manner results in maxi- 
mum support of the shield and minimizes thermal stresses. 
To be an efficient heat sink for the heating history under 
discussion (¢ max ~100 Btu/ft? sec, ¢ ~150 sec), a material 
should have a high specific heat, a high melting temperature, 
a moderate thermal conductivity and not be structurally 
brittle. A comparison of the heat sink capabilities of several 
materials for equal weight is given in Fig. 3 (prepared from 
information made available to the author by Brunner). 
The temperatures shown with the various metals (except 
beryllium at 1200 F) are several hundred degrees below their 
melting points whereas that shown with graphite corresponds 
to its sublimation point. Hence, the relative potential of 
these materials as heat sinks for the heating history under 
discussion is indicated. Beryllium and graphite are the most 
efficient, whereas nickel, molybdenum and copper are rela- 
tively inferior at these lower heat fluxes. Bervllium’ at a 
temperature of 1200 F still retains its superiority over all 
other materials except graphite (6720 F). 
Fig. 3 shows the relationship between the maximum heat 
flux (the peak re-entry flux) and the overall re-entry heating 
* Numbers in parentheses indicate References at end of paper. 
8 At present there is some concern about the structural brittle- 
ness of beryllium above 1200 F. In this discussion graphite, also 
a brittle metal, has been assumed to be intimately attached to an 
adequate structure throughout its temperature range. Beryl- 
lium, for use at temperatures when it is brittle, can be considered 
in a like manner. 
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Fig. 2 Illustrative construction—passive heat sink 
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Fig. 3 Comparison of heat sink capability of several materials 


time which heat sinks of the specified materials can accom- 
modate without the heated surface temperature exceeding 
the specified limitations. Fig. 3 applies to moderate weight 
heat sinks (~10 lb/ft). 

For the heating pulse under consideration (¢ max ~100 
Btu/ft? see, ¢ ~150 sec) the heat sinks shown in Fig. 3 are 
essentially finite slabs. Consequently, the temperature of 
the rear surface (opposite the heated surface) may assume 
intolerable levels considering the requirements of support 
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Fig. 4 Illustrative construction—ablation 


structure and internal components. Therefore insulation will 
be required behind the heat shield. The selection of the 
insulation and the placement of it relative to the heat sink 
and support structure must be done on a system optimization 
basis. Placing the support structure next to the heat sink 
(Fig. 2b) will result in greater structural weight but makes 
possible the use of low density insulation for limiting internal 
compartment temperatures. On the other hand, placing 
insulation between the heat sink and the structure (Fig. 2c) 
will result in lower structural weight but higher insulation 
weight, since the insulation will have to transmit the re-entry 
decelerative loading from the heat sink to the support 
structure. Optimization of a heat sink system as the primary 
re-entry heat protection system for a drag decelerated, 
manned ballistic vehicle, therefore, lies in two basic considera- 
tions: Selection of the best heat sink material considering the 
heating history, and tradeoffs between structural tem- 
perature and insulation requirements. 

In Fig. 4 is shown an ablation shield bonded to a low tem- 
perature structure. This type of shield serves both to remove 
some of the input heat by vaporization of the outer surface 
material and to insulate the internal structure against the 
remaining heat (Figs. 4a, b). A lighter shield construction 
would consist of a less than semi-infinite slab thickness 
ablating shield combined with an insulation layer of lower 
density than the insulating portion of the simple ablation 
shield. Such a configuration (Fig. 4c) offers substantial 
weight economies. However, since most of present low 
weight insulation materials possess low strength and since 
the insulation layer would have to support loading, suitable 
support of the insulation layer must be provided. Here also 
it is evident that for a given application an optimum weight 


can be achieved by considering various materials at appro- 


priate temperatures for the internal structure together with 
the corresponding insulation required. In the ablation type, 
however, since it is desirable to have a bond or an intimate 
support between the ablation layer and the insulation or 
main structure, it is important to select materials for the main 
structure which avoid thermal expansion incompatibilities 
between it and the ablation layer. A bond is considered 
necessary in the case of a continuous (nonsegmented) ablation 
shield because most of the efficient ablating materials are 
highly brittle and are not safe for primary load bearing func- 
tions. Table 1 (contributed by Abbate) shows some of these 
low fracture strains. At present, adhesive bonds are avail- 
able for operation up to 500 F for short times. Some tech- 
niques utilizing a combination adhesive-mechanical bond are 
being developed for bond operation up to 1000-1200 F. For 
operation at bond temperatures of 2000 F which might be 
present at the inner face of a thin graphite ablation shield, 
suitable bonds are nonexistent at the present time. 

In the absence of an adequate bonding technique for high 
bond temperatures or in the event that a continuous brittle 
shield, even though suitably supported, has excessive stresses, 
a shield composed of lateral segments (tiles) can be used to 
advantage, since the segments, if small enough and not too 
thick, can be relieved of the major thermal and load stresses. 
This also permits the use of a less rigid understructure, sinc¢ 
cracking of the shield due to deformation of the understructure 
could not occur in the segmented shield. The individual 
segments can be attached to the main structure by minimum 
restraint mechanical fastening. The presence of many breaks 
in the shield surface due to the clearance between the seg- 
ments, however, can produce turbulence and excessive heating 
in the region of the edges and corners of segments causing 
severe erosion. This deterioration of the surface contour 
with the attendant increasing turbulence and heating can be 
intolerable. This might be alleviated by making exposed 
edges blunt rather than sharp or by filling the clearance with 
a suitable material. 

As previously stated, a glide entry vehicle is exposed over 
its major area to a smaller heating rate and a much larger 
total heat input than a ballistic entry vehicle to perform the 
same mission. A heat sink is undesirable for this type of 
heating because of the large heat capacity required (and hence 
excessive weight), whereas the heating rate is generally too 
small for effective ablation-type shielding. Radiation cool- 
ing, on the other hand is ideal for use here since with the 
low heat flux the equilibrium temperatures in the presence 
of high emissivities are within the operating temperatures 
of available materials and most of the incoming heat is radi- 
ated away from the structure. A plot of the heat fluxes 
which can be radiated from a surface vs. surface temperatures 
is given in Fig. 5 for several emissivities. Indicated also are 
materials suitable for radiation surfaces. 

Shown in Fig. 6 is a cross section of a possible radiation 
cooled structure for the fuselage or cabin area. The outer- 


Table 1 Ablation shield materials +) 
= Coeff. of 
Tensile Elastic Fracture thermal 
ultimate, modulus, strain, Density, expansion, 
Material psi X 107% psi X 10-6 per cent Ib/in.* in./in. deg F X 106 
phenolic nylon 0.45 0.044 
phenolic “‘leached’’ glass 2.30 0.058 
quartz 10.0 0.07 ~0.079 9.30 
Teflon polymer* 0.10 5.00 0.080 
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surfaces would have a high and low emissivity, respectively. 
Beneath this is a layer of insulation, primarily for shielding 
against heat transferred from the inner surface of the radiat- 
ing layer. Beneath this is the main load bearing structure 
followed by a layer of insulation. Since, as is evident from 
Fig. 5, the temperature of the radiating layer is quite high 
whereas that of the inner structure is considerably lower, a 
great deal : differential thermal expansion is to be e xpec ted, 
especially in fully circular structures like internally pres- 
surized ilies. Restraining the hot outer shell against ther- 
mal expansion may cause it to buckle if it is sufficiently ductile 
or fracture if is brittle. Buckling can result in excessive 
aerodynamic losses or heating. This problem may be 
reduced to within tolerable limits by restraining the outer 
shell by means of closely spaced supports. This produces 
sniall uniform buckles having low amplitude over the surface. 
A freely expanding outer shell on the other hand may buckle 
under the aerodynamic pressures or be subject to flutter. 
The solution of this problem appears to lie in some of the 
folowing approaches: 

1 Minimize differential expansion by selecting low ex- 
pansion materials for the hotter structure and high expansion 
materials for the cooler structure, and by adjusting the heat 
transferred into the inner structure in accordance with the 
maximum tolerable differential thermal expansion (progres- 
sive stepdown of temperature). 

2 Spread out a thermal buckling tendency by attaching 
the hot shell by means of closely spaced cellular supports, 
thereby reducing aerodynamic losses and excessive heating. 

3 Allow unrestrained expansion of the hot shell only 
when it is capable of withstanding the loads applied to it. 

4 Segment the hot shell, if no other solution is adequate. 

In order to examine some of the heating problems corre- 
sponding to a radiation cooled structure in more detail, let 
us consider Fig. 7 (prepared from information made available 
to the author by Castelli). This shows the heating of the 
side walls of a vehicle experiencing the ballistic heat flux 
shown in Fig. 1, followed by a relatively long and gradual 
drag chute descent through the atmosphere (total time 1600 
sec). Two radiation wall constructions are shown, configura- 
tions A and B. These are identical except that A has an 
effective emissivity factor Fe = 0 on the interior surfaces, 
whereas B has a corresponding effective value of Fe = 0.20. 
It is seen that the outer wall in each case responds quickly to 
the heat flux. The inner wall of A rises very slowly to a peak 
of 240 F, but the inner wall of B rises rapidly to a peak of 
1000 F and only gradu ally decreases. In each case the 
temperature of the inner wall actually exceeds that of the 
outer wall for the greater portion of the descent. Thus a 
reversal of the thermal stresses of the outer and inner walls 
takes place, the outer and inner walls initially being in com- 
pression and tension, respectively, and then for the greater 
time reversed. The smaller temperature differences between 
the inner and outer walls of configuration B during the 
heating time make B appear more favorable than A with 
respect to maximum thermal stresses in the walls. On the 
other hand the temperature of the inner surface of the insula- 
tion layer (considered perfectly insulated at the surface) 
in configuration B (500 F peak) makes B inferior to A with 
respect to maintaining comfortable cabin temperatures. The 
foregoing example illustrates the necessity of minimizing the 
flow of heat into the interior of a relatively long time entry 
vehicle. It shows that for effective radiation cooling: (a) 
Very low interior surface emissivities must be maintained 
between the walls for substantially long times, (b) efficient 
radiant heat insulation is needed in the event that surface 
emissivities cannot be maintained due to oxidation, etc., or 
are not producible, and (c) conduction paths between outer 
and inner walls must be minimized. 

This example serves also to illustrate basically some of the 
heating problems of a glide vehicle such as indicated in Fig. 1. 
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Fig. 5 Surface temperature vs. radiative heat flux 
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Fig. 6 Illustrative construction—radiation cooling 
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Fig. 7 Heating of a radiation cooled structure during re-entry 
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In order to minimize thermal stresses arising in mutually 
attached inner and outer walls having widely different tem- 
peratures, a wide variety of materials must be considered. In 
doing so, thermal incompatabilities can be minimized. 
Especially in the presence of the temperature reversal effect 
demonstrated in Fig. 7, such an approach will be required. 
A selection of materials suitable for the walls of, say a radia- 
tion cooled vehicle, can be made using the information in Figs. 
8 and 9. In Fig. 8 it can be seen that for operation at the 
high outer surface temperature, materials like molybdenum 
and graphite can be used. Their low ratio of strength to 
density, however, renders them unattractive for use as 
main load bearing structures. For such use beryllium, 
inconel X and titanium among other materials appear quite 
attractive for temperatures around 1000 F. Fig. 9 shows 
that, fortunately, materials necessary for the high tempera- 
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Fig. 8 Strength to density ratio vs. temperature for several 
re-entry vehicle materials 


BERYLLIUM 
INCONEL X 
15 
MOLYBDENUM 
2 
z 107 
< 
v4 AL TITANIUM 
5 (6AL-4V) 
= 
GRAPHITE 
(RT-0003) 
7 FUZED QUARTZ 
= 
1 fe) 3000 


1000 2000 
TEMPERATURE - °F 
Fig. 9 Comparison of thermal expansions for several re-entry 
vehicle materials 


ture use have low thermal expansion. This seems charac- 
teristic of many other high temperature materials like cc- 
ramics, quartz and tantalum. On the other hand, intermediate 
temperature materials such as beryllium, 18-8 stainless, tita- 
nium and inconel X have higher thermal expansions. Finally, 
low temperature materials like aluminum and magnesiuin 
have the highest thermal expansions of all. This suggests 
that by using the right material for the hot and cooler walls 
it might be possible to limit thermal stresses to allowable 
values. In Fig. 10 are shown the thermal stresses occurring 
in several double wall constructions utilizing graphite (0.125- 
in. thickness) as the hot wall. Aluminum, titanium and 17-7 
PH steel serve as the inner wall. The thickness of the inner 
wall for each material was determined from the same bend- 
ing buckling requirement and found to be 0.029, 0.042 ani 
0.036 in., respectively. Graphite was chosen in order t» 
study the use of a brittle material under conditions of high 
thermal expansion. The graphite was assumed to be at 
constant temperature of 2240 F, and the temperature of the 
inner wall was varied. The walls were assumed to be re- 
strained by each other. For each combination a tempera 
ture is reached at which no thermal stresses exist in eith« 
wall. These are 400 F for aluminum, 780 F for 17-7 PH an 
830 F for titanium. Below these temperatures the graphit: 
is in compression, and the inner wall in tension. Abov« 
them, the situation is reversed. Considering that graphit: 
breaks at 9500 and 2600 psi in compression and tension, re 
spectively, it is seen that graphite would not be overstresse: 
provided the aluminum wall did not exceed 550 F and th: 
titanium and 17-7 PH did not exceed 1000 F. The stresse 
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Fig. 10 Thermal stresses in several double wall structures 
consisting of a hot outer wall and a cooler inner wall whose tem- 
perature varies 
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mh might be a design consideration at the higher temperatures 
ita- when they are under compression. 

lly, The development and testing of materials and necessary 
um protective coatings for operation at very high temperatures 
asts is being conducted on an accelerated basis. Manufacturing 
alls techniques of welding, brazing, diffusion bonding and flame 
ble spraying of ceramics, to mention only a few, are being de- 
ing veloped to permit fabrication of these materials for use at 
95. these required high temperatures. ig 


ner Meteoroid Impact 
ah Recent investigations of meteoroid bombardment of space 
ee vehicles have led to quantitative estimates of the structural 
igh material required to reduce the probability of penetration of 
E the vehicle wall to within conservative design limits. In 
the Fig. 11, derived from (3), is presented the estimated thicknesses 
r o! various structural materials required for reasonably safe 
ra design for a vehicle having an external surface area of 100 ft?. 
he Ii is evident that materials with a high elastic modulus are 
uncl tle most effective meteoroid shields for a given thickness. 
ite For a given weight, however, a material like beryllium will 
& be greatly superior to stainless steel and molybdenum, 
‘its because of its superior ratio of elastic modulus to density. 
re Possible wall constructions for protection against meteor- 
sa oids are shown in Fig. 12. Fig. 12a (from an idea attributed 


th: to Whipple) shows a multilayered construction in which 
the outer layer or layers, which provide the major heat 
protection, also provide the meteoroid protection. The 
inner layer is the main structural shell which also serves as 
the pressurized cabin wall. If the outer layer were punctured 
by a meteoroid, the resultant fragment should impinge on the 
inner wall over a wider area than that of the original punc- 
ture. This should make penetration of the inner wall un- 
likely. 

The apparently high efficiency of beryllium as a meteoroid 
shield suggests the use of thin sheets or foils of beryllium 
placed between the heat shield and the pressurized cabin 
wall (Fig. 12b) whenever the heat shield is inadequate for 
meteoroid protection. Such sheets could also help protect 
the inner structure from heat radiated from the heat shield. 

In any event, schemes for automatically detecting and 
sealing a penetration once it occurs should be perfected. 
This is emphasized by the fact that for a ratio of the volume 
of pressurized cabin to the area of the punctured hole equal 
to 10,000, atmospheric pressure is reduced to 3.0 psi within 
30 sec, and that the time limit of human operation at such low 
pressures is about 15 sec (6). 

In addition, the pressurized cabin shell should be con- 
structed of materials which have high resistance to catas- 
trophic bursting in the presence of a crack or hole. A com- 
parison of the bursting resistance of several structural alloys 
is given in Fig. 13 from (4). cae. 


Landing Phase 


Recovery of the vehicle may be accomplished by means of a 
drag chute or retrorocket decelerated landing on land and 
water for ballistic entry types, and more or less con- 
ventional type of landing on land or water for glide vehicles. 
The latter could also employ drag devices or retrorockets to 
assist in landing. Ballistic vehicles will experience high 
impact loadings, particularly in the case of ground impact. 
The shocks transmitted to the occupants can possibly be re- 
duced by designing a sufficient length of collapsible shell 
structure forward of the occupant mounting point. Such 
design should allow for: (a) Cushioning the occupant against 
the very short time prebuckling elastic shocks, (b) limiting 
relatively long time deceleration force applied to the human 
body during collapse of the frontal structure to within the 
human shock tolerances, and (c) maintaining the supine 
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induced in the inner walls are quite low, although buckling 
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Fig. 14 Rigid body impact decelerations of a blunt nosed vehicle 
during water entry 


body attitude and the broad distribution of forces of the 
body during the full deceleration interval. Recovery of the 
vehicle on water, on the other hand, precludes structural 
failure which would cause leakage of water into the cabin. 
Fortunately, the deceleration forces associated with moderate 
impact velocities on water can be kept within human toler- 
ances without structural collapse by proper design of the 
vehicle frontal shape. 

The landing of glide-type vehicles presents structural 
problems similar to the landing of present-day high per- 
formance fighter aircraft, except that landing velocities and 
attack angles may ‘be higher. 

We now consider briefly the recovery of a representative 
low weight to drag ballistic vehicle in water. By means of a 
drag chute or other drag devices, the flight velocity of the 
vehicle can be reduced sufficiently to permit safe landing on 
water. In Fig. 14, the rigid body decelerations of a blunt 
nosed vehicle impacting vertically on water with an initial 
velocity of 50 fps are plotted vs. the penetration depth ratio, 
d/h,. The deceleration pulse rises to a peak value of 23 g 
and falls below 5g in a total elapsed time slightly greater than 
0.1 see. In addition, the rate of increase of deceleration to the 
peak value of 23 g is less than 800 g per sec. Medical author- 
ities believe that the measure of severity of shock loading on 
human beings is the combination of the peak g and the rate 
of rise to the peak. The maximum ability to withstand shock 
loading depe ott on proper body support and orientation to 
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the loading. The best positioning of the body is supine with 
the maximum distribution of deceleration forces over the 
rear portions of the body. Shock tests with human subjects 
so supported indicate that a rate of rise of 1000 g per sec to 
a peak of 30 g is tolerable whereas a rate of 1500 g per sec 
to a peak of 40 g results in unconsciousness (5). Thus the 
decelerations shown in Fig. 14 are within the safe human 
tolerance. 

Both the peak and the rate of rise of deceleration during 
water entry are very dependent on the frontal shape of tle 
vehicle. In general, for a given vehicle mass, velocity and 
cross-sectional diameter, the greater the flatness of thie 
frontal surface the higher is the decelerative peak and rate 
of rise. Human shock tolerances, therefore, can easily he 
exceeded if sufficient attention is not given to the front.l 
shape for the water entry condition. 

The foregoing discussion does not consider the effect of 
surface waves or of a lateral velocity component such :s 
might result from lateral drift or swinging of the drag chut:- 
vehicle system. These factors, of course, must be consider¢ 
in accordance with expected overall water entry conditions 
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Mideourse Guidance 


on hem MAJOR problem of celestial navigation is the selec- 

; tion of an orbit which will intersect a given target orbit 
nd in space-time. In the case of present Earth satellites and 
for heavier Earth satellites of the near future, propulsion and 
technology limitations and/or mission requirements may 
dictate that near-circular satellite orbits be employed. 
Whether for purposes of satellite destruction or for the trans- 
fer of material to space stations, there will be a desire to 
rendezvous with these circularly orbiting satellites. 
thi A typical planned interception may be divided into four 
stages: The initial or launch guidance stage; information 
stage, i.e., determination of post-launch phase-space coordi- 
nates of satellite and intercepting rocket; calculation and 
employment of midcourse propulsive corrections; and the 
terminal guidance maneuvers. 

1 It is the function of midcourse guidance to place the inter- 
cepting vehicle in close proximity to the satellite at some 
future time so that the terminal guidance measures may be 
applied. That midcourse corrections will be needed can be 
safely predicted, not only due to launch errors, but due to 
effects of those physical space phenomena known to exist, but 
whose magnitudes have not as yet been accurately measured. 

Concerning the information stage, we note that the prob- 
lems associated with satellite tracking and the attainable 
accuracies have been discussed elsewhere (1).4 It is well 
known (2) that the determination of the satellite’s position at 
three distinct times (assuming good accuracy) is sufficient to 
vield the position and velocity of the satellite at the mean 
observation time, and hence, by Newtonian mechanics, its 
entire future path is known. As a starting point we shall 
assume that such phase-space coordinates can be determined 
by a suitable data gathering system. 

In actuality, the Newton-Kepler problem does not apply 
with complete accuracy to Earth satellite orbits. For pre- 
cision the orbital perturbations due to atmospheric drag, 
oblateness of the Earth, the sun and moon force fields, etc., 
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Midcourse Guidance Problem 


should be included. These corrections are also discussed in 
(2). Midcourse guidance as defined here, however, does 
not require inclusion of this level of precision. ; 

The particular guidance scheme which should be used for 
midcourse guidance will depend upon the type of propulsive — . 
system being employed by the rocket, i.e., whether of the — r 
sustained or impulsive type. Moreover, the appropriate — 
mathematical expansions and formulas will depend upon the ‘ 
magnitude of the anticipated corrections. We shall assume 
in this study that impulsive corrections are to be determined | 
for the interception of an Earth satellite by a rocket vehicle. 

The satellite is moving in a circular orbit, and the rocket’s | 
present orbital course is co-planar and only slightly in error. 
The orbital corrections may then be assumed small in magni- — 
tude. 


Terrestrial or Single Central Force Guidance 
Problem 


In the terrestrial problem, the equation of path is the com- 
mon conic section 
a(l — €?) 


1 + €cos (0 — A) 


(provided the energy of the vehicle is nonzero, for then the 
path is parabolic). 
The Newtonian differential equations of motion for a 
planetoid in a two-dimensional central force field require 
three constants of motion to be specified for a unique deter-_ 
mination of the path. Two are required for the radial mo- — 
tion, and only one determines the angular momentum, a 

constant of the motion for the system. 

The problem of midcourse guidance is solved when there 
are given the required changes in two orthogonal components 
of velocity such that the new trajectory will be a desired one, 
These changes in velocity are affected by a change of the 
constants of motion (determined from initial conditions), 
and since there are three such constants, two are chosen 
arbitrarily and the third is expressed in terms of these two. 
This choice nuw serves as a set of generalized coordinates 
from which the velocity changes are to be computed. 

Conversely, any given point in phase space must determine 
these orbital parameters, and in the case of the conic section | 
[1] the measured coordinates at any phase-space point in 
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transit will yield the constants of the motion via (3) 


v2 
= 1 aint y 


a = pr/(2p — rV?) 


(rV2/p) sin? — 1 
6) = 8 — cos [4] 
V1 — (rV2/u)[2 — (rV2/p)] sin? 


Here r, 6, V are the rocket’s radial and angular positions and — 


tangent orbit | 


its velocity, respectively, y is the “heading,”’ which is deter- 
mined from the ratio of transverse to radial velocities by the 
equation 

tan Y = Vo/V; 


and yw is the Earth’s gravitational constant. 

The time required for a rocket in this central force field, to 
travel from its present radial position r to the satellite’s 
radius r,, is given by 

2 Ir’ = Tm 
Vil ae ae =r 


(provided r,, is in the domain of r’). The following is a sketch 
of the derivation of Equation [5]. It is noted that 


dt = ds/V 


V is evaluated from the energy equation 
V2 
The differential arc length can be written 
ds = V1 + rX(d0/dr)? dr 


dr r2esin (0 — &) 


Eliminating angular coordinates reveals that 


~ ate — (r — a)? 


The subsequent substitution yields _ 


dt = = 
~ 


from which formula [5] follows by direct integration. Tee 


If the corresponding time for the satellite to travel an angu- 
lar distance from its present position ¢ to some @¢, is sub- 
tracted from the foregoing, a time difference is obtained 
which becomes a useful quantity when the satellite’s final 
angular position ¢» is required to coincide with the rocket’s, 


= 0... 
The angle F can be written in terms of the previous rela- 
tions as 


[6] 
=r 


(@ — ¢) is the present angle observed at the Earth, between 
the rocket and the satellite. In order for a space-time inter- 
section to occur between these two objects, the angle F must 
vanish. Thus F represents the angular difference in posi- 
tion of the satellite and rocket, that will result should the 
rocket be allowed to continue on its present orbit out to 
satellite radius. 
If F does not vanish, guidance is required, and the problem 
becomes one of determining what variation in the generalized 
coordinates a and € will result in the vanishing of F. 


From the equation of path it follows that rm +A ir 


class of all elliptical orbits \ 
with negative values of energy | 
/ | \ \~« of zero energy 
envelope starts from tangent orbit and | 
asymptotically approaches the parabola 
/ 


circular orbit 
of the satellite 


envelope spirals about the earth 
starting from the tangent orbit 


Fig. 1 Graphical illustration of the space-time restraining 
trajectories 


The variation will take place in a system with two degrees 
of freedom; hence the constraint imposed by F = 0 will be 
satisfied by a number of trajectories (Fig. 1). An additional 
constraint yielding a unique solution to the problem may be 
imposed. In particular, two alternatives based on power 
requirements are considered here. 


This alternative is accomplished by that orbit whose 
apogee is tangent to the satellite’s orbit at the point of inter- 
ception. This may be a desirable constraint because the time 
spent in the vicinity of the satellite is maximized. More- 
over, at apogee it would be possible to impart a final velocity 
thrust which would superimpose the rocket vehicle’s orbit 
on that of the satellite and the two vehicles would circle the 
Earth in permanent juxtaposition. This final thrust would 
always be needed because the rocket’s velocity at apogee is 
not equal to the satellite’s velocity. Such a maneuver may 
be preferable to nontangential interception followed by an 
attempt to change orbits, for in the latter case the velocity 
component perpendicular to the satellite’s orbit would have 
to be extinguished. 

The requirement that the trajectory’s apogee be tangent to 
the satellite’s orbit may be expressed by the following equa- 
tion 


Tangential Interception 


G=rn—ail+e)=0 [7] 


Presumably G will not be zero, and it is a measure of how 
far the turning point P of the orbit has missed tangency with 
the satellite. In theory then, Equations [6 and 7] can be 
solved for a and ¢, requiring F and G to be zero. With these 
values of a and ¢ the required velocity components can be 
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Now, under the assumption of small corrections 


FP 
‘4 + da, € + de) Fla, 6) + da + de = 0 


similarly for G(a + da, € + de). We compute the correc- 
tions da, de which will make F and G zero from the equations 


Solving 


(OF /d€)G — (dG/d€)F 


~ (OF /da)(0G/de) — (OF/d€)(dG/da) 


da 


(0G/da)F — (OF /da)G . 


de = (SF /0a)(0G/de) — (OF /de)(0G/oa) 


These formulas are true in general. The partial derivatives 
are easily computed for the F and G given, and they are 
collected at the end of this section on the terrestrial or single 
central force guidance problem. (They have not been sub- 
stituted here, for they appear to reveal no new physical 
consequences.) 
The corresponding incremental velocity changes are 7 


_# 


Minimum Midcourse Guidance Power 


This alternative is accomplished by choosing the impact- 
ing orbit requiring the minimum expenditure of fuel during 
midcourse correction. 

A similar minimization problem is well known in the liter- 
ature. This problem has been formulated in the following 
manner. 

Given the phase-space coordinates of two points in a 
gravitational field, and the equation for the burning rate of 
the fuel, what course is to be followed between these two 
points such that minimum fuel is expended? Mathemati- 
cally, this is a problem in the calculus of variations in which 
two points in phase-space are given, and the relative mini- 
mization of the integral for fuel consumption over all co- 
terminal curves is considered. This solution merely will 
give the path in which a minimum amount of work is done 
against the fields, with most of the work going into changing 
the kinetic energy of the rocket. This important problem 
is distinct from the problem that is being considered here; 
namely, if the rocket is off course, what revised orbit shall be 
taken at that point which will evolve successfully to the 
satellite intercept and how shall the alterations be brought 
about with a minimum expenditure of fuel? 

This minimal condition must be a scalar, for ultimately it 
is the change in mass of the fuel which is to be minimized. 
The potential energy of the fuel, be it chemical or nuclear, 
is proportional to its mass, and any change in the mass of the 
fuel must reflect itself in an equivalent change in kinetic 
energy of the rocket. The momentum given to the field can 
be neglected for impulses of small time durations. 


A (lode | oa oe da 

oa 
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The change in kinetic energy of the rocket is given by 


dK = 2V-dV + dV-dV 


First-order terms are not sufficient as a measure of this 


change, for a maneuver consisting of 
VodVo = -V, dV, 


would cause the first-order term to vanish; however, this 
would not be a measure of the power used in the transition. 
On the other hand, the complete expression cannot be the 
correct one either, for in the case of a mere rotation of the 
velocity vector, power is expended, but the energy of the 
initial and final states is the same; moreover, the complete 
expression would be a measure which is velocity dependent. 
An attractive scalar to minimize is one that measures only 
the power expended in the transition, and is not proportional 


tothe velocity 


dk = dV-aV. 


This quantity is to be minimized subject to the constraint of 
space-time intersection, i.e. 


F(a, €, da, de) = 0 


Oa O€ 


E Ve vel 


Consistent with the previous assumption, Taylor’s theorem 
is used to the first order to measure F 


oF oF 
F(a, €, da, de) = F(a, €) + — da + — oe 
oa de 


Here a and ¢ are given parameters computed from the present 
phase-space coordinates, and the variable quantities da and 
de are to be chosen on the basis of the two foregoing condi- 
tions. 

By the method of Lagrange multipliers m must be elimi- 
nated from the three equations (following), and then the in-— 
dependent variables da and de must be obtained. 
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The solutions are 
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aaa Angles are measured in degrees. 
a Distances are in satellite radius units. 


Fig. 2 Required eccentricities for trajectories which contact the 
satellite at the apogee 


semi-mojor axis 


6 8 
eccentricity (€) 
Fig. 3 Navigation plot of the guidance parameters for the range 
r/tm = 0.5 


1.0 


Region 1—maximum range is too short, and range never assumes 

r/rm = 0.5; region 2—maximum range is sufficient, but range never 

assumes r/rm = 0.5; region 3—maximum range is too short, but 

range does assume r/rm = 0.5; region 4—maximum range is suffi- 
cient, and range does assume 7r/rm = 0.5 


oVe\? 
a= + -2 

| 


Oa 

The differential changes da and de vanish when F vanishes. 
Physically this is not only to be expected but lends support 
to postulate that dV-dV is an appropriate quantity to he 
minimized. If the complete change in kinetic energy weve 
taken as the minimal condition, it would contain terms line:r 
in da and de, and thus nonzero values of da and de wou'd 
appear when F is zero, i.e., even when the principal con- 
straint of space-time intersection is met a maneuver wou'd 
be called for. This is an absurdity since only the nonvanis).- 
ing of F should dictate the execution of the maneuver. 

The partial derivatives required for explicit determinaticn 
of the velocity changes are 


Vu Late 1 — [(r — a)/ae]? — [(r — a)/ae)® 
atl 
aev/ 1 — [(r — a)/ae]? 
a3 
Vi 
= —(e+ 1 


r — a(l — e€)? | 
— [(r — a)/ae}? 


2V,a? 


E = = 
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Numerical Example of Satellite Interception 


This section displays sample impacting trajectories and a 
semi-analytic, or graphical solution to the midcourse guid- 
ance problem as applied to the space vehicle in a single force 
field. It is assumed here that the satellite orbits the Earth 
in a circular path. If the vehicle is to intersect the satellite 
in space-time, it has been shown to be necessary for a and € to 
satisfy the relation 


(2) 


This expression is obtained from 6 when all distances are 
normalized by the satellite’s radius r,, and the satellite's 
angular velocity V u/tm is taken as the unit of angular 
velocity. 

A convenient coordinate system here is one with an origin 
at the Earth’s center and the z axis through the satellite, 
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Fig. 4 Velocity component contours as a function of focal angle 
for satellite interception when r/rm = 0.5 


‘so that the system rotates with the angular velocity of the 
satellite. The polar coordinates of a particle in this system 
are then (@ — 0,1:/rm). Trajectories of different eccentrici- 
ties which contact the satellite at their apogees are shown in 
this coordinate space in Fig. 2. The trajectories are relative 
to the Earth satellite axis and depict how the rocket positions 
would appear to an observer on Earth who faced the satellite 
at all times. 

The set of conditions 


describe a region in (a, €) space within which these trajectory 
parameters must lie in order to intercept the satellite. This 
region is a function of the distance to the vehicle and is 
shown in Fig. 3 for r/rn = 0.5. Within this region the 
(a, €) combinations of practical interest are those in the 
shaded subregion which correspond to ellipses of high eccen- 
tricity ; those not wasteful of power. 

Now at any point (@ — 6, r/rm) in space, all values (a, €) 
obtained from Equation [8] and satisfying conditions [9] 
will allow space-time intersection with the satellite. Since 
the variables under control of the vehicle are the velocity 
components, it is desirable to transfer the consideration to 
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Oscilliscope showing the 
steering dot and velocity contour 


Range Focal Angle 
Indicator Indicator 


- @ 


Fig.5 Instrument layout for presenting the guidance information 


where JV, is satellite velocity. 

The dashed envelope in Fig. 4 represents the inner region 
of Fig. 3 as plotted in V,, Ve space, at r/rm = 0.5. 

For any angle (9 — ¢) measured at this range, the terminal 
point of the velocity vector must lie upon the appropriate 
contour if satellite interception is to be accomplished. If it 
does not, sufficient rocket power must be applied to bring 
the velocity vector to the required contour. 
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in Fig. 5. Here the single velocity contour, ABC (as a func- 
tion of the measured range and angle) is continually gener- 
ated and can be presented on the face of a scope. The 
measured velocity is also presented as a steering dot. If the 
missile is on a proper trajectory, the dot should remain on the 
contour throughout the flight. If it is off, the dot must be 
“steered” onto the contour. Any point on the contour will 
provide an interception; two points are of particular interest: 

1 Point A: This point provides a tangent interception. 

2 Point B: This point, nearest to the dot, requires the 
least midcourse guidance power. 

It is of further interest to note that the steering dot ini- 
tially does not have to fall within the outlined region. A 
point to the left of the region merely means that the rocket 
will not reach satellite’s radius without the application of 
additional power. 


Nomenclature 

im = radial distance between the Earth and the satellite 

r = present radial position of the rocket 

6 = present angular position of the rocket 

Om = angular position of the rocket when its radial position 
is Tm 

¢ = present angular position of the satellite 


N ALL guidance systems it is necessary to maintain knowl- 
edge of missile angular orientation. In a command guid- 
ance system, maintaining this knowledge at the command 
point poses a special problem because of the inability of track- 
ing devices to measure orientation of a remote missile (1).? 
Some degree of cooperation from the missile is necessary, and 
this may be achieved in a variety of ways. One method, 
which is sometimes necessary but substantially increases 
missile cost and complexity, makes use of a complete missile- 
borne inertial reference capable of transforming received 
commands phrased in inertial coordinates into missile pitch 
Received Feb. 12, 1959. 
1 Engineering Specialist. 
? Numbers in parentheses indicate References at end of paper. 
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A further example of the guidance technique is displayed 


om = angular position of the satellite when the radial posi- 
tion of the rocket has reached rm 

w = angular velocity of the satellite when circular motion 
is assumed 

a, ¢, % = semimajor axis, eccentricity, and angular orientation 
of the elliptic conic section 

F = focal angle between the rocket and the satellite 


measured at the Earth, when the radial position 
of the rocket has reached the value rm 

G = distance from the turning point P of the orbit t« 
the trajectory of the satellite measured from thi 
Earth along a common radius 


if = turning point or maximal point of the orbit from the 
Earth 

V; = radial component of the velocity of the rocket 

Vo = transverse component of the velocity of the rocket 

V = speed of the rocket 

m = gravitational constant of the Earth 

y = are tangent of the quotient of V; to Vo; see Fig. 6 
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Some types of command guidance systems utilize tracker data plus the missile roll autopilot con- 
trol equation to compute missile angular orientation at the command point. 
tracker errors, autopilot transient errors and angles-of-attack result in incorrect command trans- 
mission and a consequent performance degradation. 
performance degradation resulting from tracker errors and angles-of-attack in two types of com- 
mand guidance systems employing the roll autopilot control equation for orientation computations. 
It is shown that much smaller command errors are expected to occur in a system in which the roll 
autopilot keeps the time integral of roll rate equal to zero, than in a system in which the roll auto- 
pilot maintains perpendicularity between one wing plane and a fixed inertial reference plane. 


In these systems 


An analysis is made of command errors and 


and yaw commands. Another method, which lessens the bur- 
den on missile-borne components and which in some cases 
functions satisfactorily, employs the properties of the roll 
autopilot control equation to deduce (at the command 
point) missile orientation. Commands transmitted to the 
missile can then be phrased directly in pitch and yaw terms. 
Unfortunately, missile orientation cannot be deduced with 
perfect accuracy by this method and a degree of guidance 
degradation is introduced. 

Command errors and guidance degradation in command 
guidance systems employing two types of roll autopilot con- 
trol equations are analyzed below. These two types of auto- 
pilot control are designated bank stabilization and roll 
stabilization, respectively. The bank stabilization roll 
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restricted to a cruciform missile configuration. 


|] in the roll stabilized system 


d= ig sin dy 


autopilot is subject to transient errors. 


tracker inaccuracy than the bank stabilized system. 


Coordinate Geometry 


line 


N = #) sec 6 


is as shown. 


ertial, missile body or Euler axis components. That is 


= UoXo + + (inertial) 
= + + Usd (body) 


| 


= + UWeN + UseZo 


(Euler) 
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autopilot maintains perpendicularity between one wing plane 
and an inertia-fixed reference plane. (In some applications 
this is called the belly-down or one-wing-vertical system.) 
The roll stabilization roll autopilot keeps the time integral 
of roll rate equal to zero. Typical equipment configurations 
(2,3) for both types of attitude control are shown schemati- 
cally in Fig. 1. For convenience, consideration has been 


The manner in which missile orientation is deduced at the 
command point is common to both bank and roll stabilization 
systems. In essence, missile orientation is completely de- 
fined by three angles, y, @ and ¢, or azimuth, elevation and 
bank respectively, regardless of whether these appear ex- 
plicitly in guidance computations. Further explanatory 
material can be found in the following section. Approxima- 
tions to azimuth and elevation of the missile longitudinal 
axis are obtained from tracker measurements by assuming 
that the missile longitudinal axis is colinear with the flight 
path tangent. Bank angle @¢ is either assumed to be zero 
in the bank stabilized system or is computed from Equation 


The direction cosines of the missile pitch and yaw axes with 
respect to an inertia-fixed coordinate system are trigonome- 
trie functions of y, @and ¢. Errors in command point com- 
putation of these direction cosines result in incorrect pitch 
and yaw command transmission with a consequent guidance 
degradation. Errors appear initially in y and @, because 
angle-of-attack invalidates the assumption of colinearity 
of missile longitudinal axis and flight path tangent; further- 
more, tracker measurements are not perfect and the roll 


It is shown in the following sections that the effects of 
angle-of-attack and tracker inaccuracy in the bank stabilized 
system differ materially from those in the roll stabilized 
system. (Angle of attack, or sideslip, and an error in meas- 
urement of direction of the flight path tangent have identical 
effects upon command errors.) Specifically, the roll sta- 
bilized system is much less affected by angle-of-attack and 


With reference to Fig. 2, the set of mutually orthogonal 
unit vectors Yo, Yo, Zo, is fixed in inertial space. The rela- 
tive orientation of a second mutually orthogonal set of unit 
vectors ¥, 9, 2 fixed to the missile body, is given in terms of 
the Euler angles y, 6 and ¢ (azimuth, elevation and bank 
angles respectively). A seventh vector N is defined by Equa- 
tion [2], and is therefore a unit vector colinear with the node 


The directions of £, N and Zp are frequently called Euler axes 
and correspond to the rotation axes of a three-axis gimbal sys- 
tem. In Fig. 2,an Xo — Yo inertial reference plane is shown as 
a horizontal plane; however, all following statements and deri- 
vations apply irrespective of %5 — Yo plane orientation, pro- 

viding the relative orientation between £, 7, 2 and Xo, Yo, Zo 


Any vector # can be expressed as the vector sum of in- 


It follows that the instantaneous missile angular velocity can 
be expressed in either of two ways as the vector sum of angu- 
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Fig. 1 Bank and roll stabilized command guidance systems 


TARGET TRACKER 
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Fig. 2 Coordinate geometry. As shown y and ¢are positive, 6 is 


negative 
larrates. That is 


@ = + ON + = [4] 


where @ is the instantaneous vector angular velocity and p, 
q, r are the missile angular velocity components along %, 7 
and 2 respectively. Specifically, p is roll rate and it can be 
shown that 


p = @ — (sin Ay 


In the roll stabilized case 


so that, in this case 
de = sin 6 dy {6] 


Command Errors 


In both the bank and roll stabilized systems, approxima- 
tions to y and @ are computed from tracker data by the 
assumption that % is colinear with the flight path tangent, 
an assumption which is generally not true because of angle- 
of-attack (or sideslip). In the bank stabilized system, 
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¢ = 0 by virtue of autopilot control, and in the roll stabilized 
system, ¢ is computed from Equation [11]. Let wy’, 0’ and 
¢’ be the apparent, or approximate values of the true Euler 
angles y, @and ¢. Then, with reference to Equation [3], if 


i= (uXo + + usoZo) [7] 


is the desired command vector at any instant of time (for 
example, % may be the desired acceleration vector), the 
ground computed pitch and yaw commands are : 


(pitch command) = 


(yaw command) = u's 

where, from the usual coordinate transformation equations (4) 
(sin sin + cos W’ sin 0’ cos + 
(—cos sin + sin W’ sin 0’ cos b’)um + 

(cos 8’ cos )uso 
= (—sin cos + cos sin sin d’)wo + 
(cos cos + sin sin 6’ sin + 

(cos 6’ sin ’)uszo 


Us = 


(Note that the normal force produced by a pitch angle-of- 
attack, rotation about the 7 axis, is directed along the 2 
axis. It is assumed that means for producing command 


accelerations along the # axis are not available.) 


ED 


\ DESIRED COMMAND 
MAGNITUDE ERROR 


Command geometry missile viewed head-on 


«© 


APPARENT DIRECTION OF 
FLIGHT PATH 


® TANGENT TO 


ANGLE 7 IS CUE TO COMBINED 
ANGLE-OF-ATTACK AND 
SIDESLIP 


These commands are executed by the missile along the 
true § and 2 directions, so that, whereas the desired response is 
(und + Us2) 


the actual response is 


+ 
If the missile is viewed head-on, looking along the longi- 
tudinal, or = axis, the command geometry will appear as 
shown in Fig. 3. In this way, the command error can be 


visualized as a rotation € and a fractional magnitude change 
6 where 


Usb 
U x» 
1, 1/2 
+ Usp? 


A preliminary comparison of the bank stabilized and roll 
stabilized systems is presented below in terms of the errors 
€ and 6 under a limited number of conditions. It is assumed 
that at some time ¢ all errors are zero and Yy = ¢ = 0. In- 
mediately after time t, the missile trims a particular combina- 
tion of angle-of-attack and sideslip designated 7, Fig. 4, such 
that the radar error is 


In the presence of this radar error, a command # = kY, is 
transmitted to the missile. Since 7 as described produces a 
change in y only, the relation between y and 7 is readily 
determined. With reference to Fig. 4 


—sin = 


Thus, from Equations [6, 8 and 10]: 
In both systems 


sec 6 sin n/2 [10] 


tho = Us = O 
uo = k 
In the bank stabilized system 
=¢'=0 
«= Kk sin sin 0 
u» = k eos 
0 
= k 
In the roll stabilized system 
ux = {—cos Wsin (sin 0) + sin sin 6 cos (W sin 6)}k 


1 O sin sin 0)}k 


The resulting error € computed from equations [9 through 


13] is plotted as a function of 9 with @ as a parameter in 
Fig. 5. These results show a rather striking superiority of 
the roll stabilized system over the bank stabilized system 
under the specific conditions chosen. It should be noted 

- that the choice ¥’ = Y(t) = 0 is not at all restrictive, since 

the direction of Yo in the “horizontal’’ reference plane is com- 

pletely arbitrary. 
both systems and does not exceed 0.15 for 7 < 30 deg. 


The magnitude error 6 is comparable in 
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Guidance Degradation Due to Command Errors 


¢ 


An intuitive appreciation for the effect of a rotation error 
e can be obtained by analysis of a simplified sampled data 
guidance model in which missile position measurements are 
made at regular intervals and a lateral translation of the 
missile is effected without delay in response to command 
immediately after a measurement is made. In the case 
considered, the missile is initially flying a straight line path 
parallel to, but displaced from, the desired straight line tra- 
jectory. It is assumed that differences in displacement of 
the missile from the desired path, measured normal to the 
missile longitudinal axis, due to changes in missile heading 
are small enough to be neglected. Heading changes of 
+30 deg are probably acceptable in this context. 

The error geometry viewed in a missile head-on perspec- 
tive is as shown in Fig. 6. A coordinate system fixed to the 
desired trajectory is oriented so that the initial displacement 
is along the abscissa. The missile vector displacement at the 
ith data sampling time is F,, and the vector command dis- 
] placement executed by the missile immediately thereafter is 
C;. Ideally C; would equal however, due to C; is 
unintentionally rotated from —£; and because of this rota- ae 
tion it is desirable to modify the magnitude of C; by a factor 
A. Thus 


8=30° 


(DEGREES) 


€ 


ROLL STABILIZED SYSTEM 


BANK STABIL:ZED SYSTEM 


— 


= 


3 0 20 30 


and if en 
Fig. 5 Rotation error « vs. 7 


then, in matrix notation 


tom tas 


a; cos € —sin €\/a; 
bi sin e cose /\b; v, AND V5 ARE UNIT VECTORS 


which can be reduced to mf i 
= M {17] DESIRED TRAJECTORY 
b; 


where 2 
1 — A cose —A sin E; = = 
A sin € 1 — A cose Cc, .s 
€ Ci-2 
Consequently 
MISSILE 
‘)=mi( 19 


where 
Fig. 6 Error geometry 


Aas 


a 
( ) = initial vector displacement 
A 


Let 


A quasi system rise time 7 can be defined as the number of 


p = (1 — 24 cose + A*)” sample periods required to reduce the ratio |F;|/|Eo| to 


0.1. It follows that 
which is the product of a scalar p and a magnitude preserving It is evident that the sequence £; will not converge if p > 1, 
rotation matrix. Therefore and further, that 7 will be minimum when p is minimum. 
Minimization of p is accomplished by adjustment of 4A, 
= pi eee [22] Equation [21]. A tabulation of 7 vs. € for two values of A, 
minimizing p for e = 0 and e = 30 deg respectively is given 
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Table 1 Rise time 7 vs. command rotation error « 


——— (sample periods) ——— 
p minimized at 
« = 30 deg 


p minimized at 
e = 0 deg 


gs 


The value of 7 listed is the integer nearest the 
[23]. 


in Table 1. 
value computed from Equation 


Conclusions 


From the preceding analyses, it is concluded that in many 
cases? much smaller command errors occur in a roll stabilized 
system than in a bank stabilized system. It has been shown 


3 The author has been unable to conceive of any case in which 
the reverse is true. 


An Approximation to the Re-Entry 
Trajectory 


that situations can readily be synthesized in which errors 
in the two systems differ by an order of magnitude. 

For example, with reference to Fig. 5, if @ = 60 deg and if 
the combined angle-of-attack and angle-of-sideslip 7 is 20 
deg, then the rotation errors are approximately 37 and 2 deg 
in the bank and roll stabilized systems respectively. Further 
reference to Table 1 shows that in this case the bank system 
rise time, 7, is more than three times as large as the roll sys- 
tem rise time. In both systems a reduction of command 
errors can be accomplished by guidance computer estimation 
of angle-of-attack from missile lift characteristics and com- 
mand history; however, this additional correction is probably 
unnecessary in a roll stabilized system. 

In addition to the roll stabilized system advantages shown 
in this analysis, it should be noted, with reference to Fig. 1, 
that the roll system rate integrating gyro (or rate gyro plus 
integrator) will never tumble and is inherently a more 
rugged device than the bank system free gyro. 
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An approximation to the re-entry trajectory of a rocket or satellite is obtained by writing an 
integral equation for the trajectory angle and using simply the re-entry angle as a trial solution. 
The integrals can be evaluated if a constant drag coefficient and exponential atmospheric density 


are assumed. 


HEN a rocket or satellite enters the Earth’s atmosphere, 

its trajectory is strongly influenced by air resistance. 
An accurate determination of the trajectory usually requires 
a numerical integration of the equations of motion. However, 
an approximate analytic solution can sometimes be useful for 
obtaining the important features of the descent. For ex- 
ample, the Allen-Eggers approximation,? which takes into 
account the effect of air resistance but neglects the gravita- 
tional force, is often adequate for calculating the maximum 
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trajectory angle does not change greatly. 
approximation gives the qualitative behavior of the descent correctly. 
weight-to-drag ratios and initial conditions are compared with exact solutions. 


The approximation is good if the re-entry angle is greater than 5 deg and if the 


However, even for trajectories which bend greatly, the 
J 
Results for a wide range of 


deceleration of the vehicle. The purpose of this paper is to 
give an improved approximation which includes the effects of 
both air resistance and gravity on the re-entry trajectory of a 
nonlifting, freely falling body. sr? 

The force due to air drag is usually expressed as 


Fp = (1/2)p(r)v?C pA 


Equations of Motion 


where 


p(r) = density of air at a distance r from the Earth’s center 
-.- (assuming a spherically symmetric, nonrotating 


atmosphere) 
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v= the speed of the falling body 
Cp = drag coefficient 
A = cross-sectional area 


The gravitational force may be written 


= mg(r.*/r*) 


where 
- m = mass of the falling body 
~ g = gravitational acceleration at the surface of the Earth — 
= Earth’s radius 


With these forces the equations of motion are 


d*r dp Fp. 
a eae 6=0 [1] 


ld 


@ = polar angle of the vector r 
i) = angle measured downward from the local horizontal to 
the velocity vector as shown in Fig. 1 


The relations 


dr/dt = —v sin 


r(d@/dt) = v cos 
may be used to eliminate the time and transform Equations 
{1 and 2] to the following equations? in v and @ as functions 


ll 
Ww 


of r 
f(r) 2} du f(r) 
dr? E 0 dr r sin 3] 
gre 
Ar) = — 4 
r%d/dr) In (r cos @) [4] 
where aw 
u = 1/cos?6 


If Equation [3] could be solved for @ as a function of 7, then 
the speed of the body could be found immediately from Equa- 
tion [4]. Unfortunately, the solution to this kind of nonlinear 
equation is not known. However, something can be learned 
about it by transforming to an integral equation in the follow- 
ing way. Let 

du 2u 


F(r)=—- - — 


dr r 
Then Equation [3] may be written 
~ 


sin sin 0(r) 


Assuming, for the moment, that sin @(r) is known as a function 
of r, we integrate this equation to obtain 


du Qu du 2uo r f(r) 
dr or - (4, exp [- sin 0(r) ar | 


where the subscript zero refers to initial conditions at some 
distance rg where the body may be considered to enter the 
atmosphere. 


8 These equations were established in an unpublished memoran- 
dum by A. D. Wheelon of Space Technology Laboratories. 
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Fig. 1 Re-entry trajectory 


Vith the help of Equation [4] it is easy to show ‘tha 


2u0 _ gr? 


cos? 2 Ay 
Integrating paar [5] we obtain 


dr 

cos? A(r) cos? r cos? Or) 

J. | - sin O(r | va | 
‘This integral form shows us how the trajectory angle depends 
on the initial conditions and on the drag through the quantity 
f(r)/sin 6. From Equation [4], we find that the speed of the 
body is given by y 


To? cos? f(r) 
= vo* exp 8 
a 
Approximate Trajectory 


Equations [7 and 8] suggest an approximation for trajec- 
tories with values of 4 which are not too small. Let us sup- 
pose that the trajectory angle # changes slowly so that it may 
be replaced by @ wherever it appears under an integral sign. 
This should be a good approximation for the initial portion of 
the trajectory. Let us also approximate the density of the 
atmosphere by 


p(r) = 
where 
kh =f—f, 
p, = density at the surface of the Earth 


Then f(r) may be expressed as 
f(r) = 
pg(CpA/W). 


With these approximations 


where A 


F 
D 
r 
| 
q 
a"? 
- 
= 
x 
: : 


Equations [7 and 8] become 


+2In + [Ei(z) — 
0 
» 
wh) = 10] 
cos 6 r 
where 
—ah 
asin 
—ahg 
asin 
Ei(z) = dz 
“ee 


Although this approximation should be valid only over the 
initial portion of the trajectory, it has been used for the entire 
trajectory in Table 1 to show how far it may be carried. The 
initial conditions are typical re-entry conditions for rockets 
and satellites and the weight-to-drag ratios W/CpA cover a 
reasonable range. The value of ho was taken to be 250,000 ft 
and the density of the atmosphere was approximated by using 

= 1/24,000 ft“! and p, = 2.37 X 10~* slug/ft®. For com- 
parison the exact trajectories, obtained by numerical integra- 
tion‘ of Equations [1 and 2], are also given. The program for 
numerical integration usually employs a variable Cp and a 
standard atmosphere, but in this case a constant Cp and an 
exponential atmospheric density with a = 1/24,000 ft~! and 
p, = 2.37 X slug /ft® were used. 

The approximation can be improved by using expressions 
[9 and 10] down to an altitude h,; where @(h,) is found to differ 
considerably from the initial angle %. Then the values », = 
v(h;) and 6, = @(h,;) can be used as initial conditions for 
calculating the remaining portion of the trajectory by Equa- 
tions [9 and 10] with 4, vp and 2 now replaced by @,, v; and 2. 
If desired new values of a and Cp may be used to give a more 
realistic trajectory. Several such steps should greatly improve 
the results. However, if an accurate trajectory is needed and 
a high speed computer is available, it is far faster to integrate 
numerically Equations [1 and 2] directly than to use many 
steps in this approximation procedure. 
entry trajectory is 


® re dr 
R= rf, do f cot 6 - 
0 To r 


By approximation [9] this is a 
re r 
R= —r cot 


2 dr 
“(Bile = {11] 
tha 
Ap = 2ge-*e/avo? sin? al, 


If the trajectory is short range, so the Earth’s curvature is 
unimportant, these approximations take a simpler form suita- 
bie to a “flat Earth.”” By letting 79, 7, and r approach infinity, 
we obtain 


where 


cos 4 
os 


‘The numerical integration was done on the IBM 704 by a 
poe used for trajectory studies i the Applied Mevriaiins 
epartment of Space Technology Laboratories. 


Pry 


The range of the re- 


The range becomes 
R = —cot % + Ao[Ei(z) — Ei(z)]}~*”? dh [14] 


or 
R= oot (1 + [15] 


a V Bo 


where 
4 
By = 1 — AoEi(z) 


ze = A/asin 4% 


As a special case we set g = 0 to remove the effect of 
gravity. Equations [12 and 13] then yield the Allen-Eggers 
approximation 


Oh) = wh) = 


The body moves in a straight line with its speed decreasing 
due to air resistance. This approximation is good for bodies 
with high weight-to-drag ratios, since such bodies are not re- 
tarded enough for gravity to produce an important bending 
of their trajectories. 

With the gravity term retained, Equation [15] for the range 
may be evaluated approximately by dividing the range of in- 
tegration into several parts. From z = % to z = 1, CpHi(z) is 
usually small compared to 1. Therefore, the quantity [1 + 
C,Ei(z)]~'/* can be expanded and the first two terms easily 
integrated. From z = 1 to z = 10, Ei(z) can be approxi- 
mated by a few quadratic sections, again making an integra- 
tion possible. Beyond z = 10, C,#i(z) is usually large com- 
pared to 1. If we neglect 1 and use the first term e?/z of the 
asymptotic expansion of Hi(z), the integral from z = 10 to 
z = 2z, (if z, extends into this range) becomes an error function. 
If we use this method to estimate, for example, the re-entry 
range of the first trajectory in Table 1 (W/CpA = 10 lb/ft?, 
% = 45 deg, vu = 8000 fps), we obtain 1.34 X 10° ft com- 
pared to the exact range of 1.49 X 10° ft. (If the body fell 
under the force of gravity alone, the re-entry trajectory would 
have been approximately a parabola of range 2.24 X 10° ft.) 
The range on a round Earth, as given by Equation [11], can 
be approximated in a similar manner. If the term (2/sin? 4) 
In (r/ro) is small compared to the other terms, the quantity 
in the square root may be expanded suitably and the logarith- 
mic term integrated to give a correction to the flat Earth 
range [15]. For example, the range of the second trajectory 
in Table 1 (W/CpA = 10 lb/ft?, 0 = 20 deg, v) = 24,000 fps) 
is 3.86 X 10° in this approximation compared to 3.98 X 10° ft 
for the exact range.® 


(2-40) /2 [16] 


Conclusions 


This approximation to the re-entry trajectory gives an ex- 
tension and improvement of the Allen-Eggers result (Eq. 
[16]). In particular, it gives a better value for the speed of 
the vehicle, it approximates the curvature of the trajectory 


and provides a reasonable estimate of the re-entry range. 


The examples in Table 1 show that, as expected, the ap- 
proximation is good whenever the changes in the trajectory 
angle are small. Even when the angle changes by a large 


5 It is difficult to say how much the exact (or model) trajectories 
given here differ from actual trajectories, since the latter are deter- 
mined partly by the dependence of drag coefficient on Mach num- 
ber for the particular re-entry body. To obtain some measure 
of the discrepancies one should expect, a number of actual trajec- 
tories were calculated by a numerical integration using the density 
of the ARDC model atmosphere (1956) and three different curves 
of drag coefficient vs. Mach number. It was found that the re- 
entry ranges of these trajectories differed by —1 to +3 per cent 
from the ranges of the model trajectories. The final speeds were 
in all cases of the same order as those of the model trajectories, 
varying from 50 to 130 per cent of those given in Table 1 in the 
columns labeled ‘‘exact.’’ 


ARS JourRNAL 


= 
> = > 
» - 
- 
- 
Ne 
> 
> 
~ 


4) 


W h, ft X 6, vX10-% 


CpA \ ft? 103 deg fps 

10 250 45.00 8.00 45.00 8.00 20.00 
200 46.25 7.98 46.24 7.98 20.07 

150 47.64 6.59 47.62 6.62 20.37 

100 55.40 1.35 54.43 1.54 53.15 

50 90.00 0.222 89.99 0.266 90.00 

0.000 90.00 0.000 90.00 0.092 90.00 

250 45.00 8.00 45.00 8.00 20.00 

200 46.22 8.17 46.22 8.17 20.06 

150 47.35 8.19 2.13 

100 48.61 7.05 48.58 7.09 20.38 

50 53.44 1.91 52.86 2.15 36.75 

0.000 90.00 0.249 89.94 0.299 90.00 

250 45.00 8.00 45.00 8.00 20.00 

200 46.22 8.19 46.21 8.19 20.06 

150 47.34 8.37 47.23 $.37 20.12 

100 48.40 8.41 48.37 8.41 20.19 

50 49.51 7.49 49.46 7.53 20.38 

0.000 2.43 2.59 62.33 2.387 28.10 


Table 1 Approximate trajectories computed from Equations [9 and 10] compared to exact trajectories obtained by numerical in- 
tegration of Equations [1 and 2] 


——Approximate— ——Exact——. —Approximate— ——Exact——. —Approximate————Exact—— 


24.00 20.00 24.00 5.00 30.00 5.00 30.00 
22.75 20.07 22.73 4.78 24.10 4.79 23.88 
14.53 20.36 14.60 10.43 4.15 9.19 5.90 
0.601 41.14 1.02 90.00 0.220 59.22 0.856 
0.153 89.99 0.266 90.00 0.000 90.00 0.265 
0.000 90.00 0.092 90.00 0.000 90.00 0.092 


24.00 20.00 24.00 5.00 30.00 5.00 30.00 
23.93 20.05 23.93 4.60 29.40 4.60 29.36 
22.94 20.13 22.94 4.37 24.66 4.40 24.06 
16.01 20.38 16.03 8.19 5.98 8.06 7.00 
1.01 32.16 1.46 89.99 0.247 51.63 0.987 
0.171 89.95 0.299 90.00 0.000 89.97 0.299 
24.00 20.00 24.00 5.00 30.00 5.00 30.00 
24.05 20.06 24.05 4.58 29.99 4.58 29.99 
24.01 20.11 24.01 4.13 29.51 4.14 29.42 
23 . 22 20.18 23.25 3.82 25.64 3.90 24.59 
17.40 20.37 17.45 6.49 8.21 7.13 8.05 
1.78 26.62 2.21 89.83 0.278 58.24 0.879 


amount (for example, from 5 to 90 deg), it gives qualitatively 
the rapid bending of the trajectory near the bottom of the 
atmosphere. The results may be improved by dividing the 
trajectory into several sections in which different values of 
ho, % and @ are used. The appropriate division can be esti- 
mated from the rate at which the angle @ is found to change. 
This approximation cannot, of course, replace machine com- 
putations when accurate results are needed. It serves 


mainly as an analytical structure for understanding re-entry 
trajectories and predicting their general behavior. 


Acknowledgment 


The author wishes to thank Dr. A. D. Wheelon for sug- 
gesting this problem. The development of this approxima- 
tion is due largely to his encouragement. 


Unstable Burning—Theoretical 
Unstable Burning—Experimental 
Grain Design 


ARS SOLID PROPELLANTS CONFERENCE 
January 28-29, 1960 


Princeton University, Princeton, N. J. 


Combustion, Gas Kinetics and Ignition 
Interrelations of Nozzle Geometry 
and Nozzle Performance -#e 


| 
of 
o 
5 
S 
— 
e 
a 
- 
= 
. 
- 
+ 
~ 


_ by the average power generated. 


energy can be stored in batteries 


pu PROBLEM of making a transfer between two orbits 
has been considered by many authors. For example, 
Hohmann (2)? showed that for a transfer between two circular 
orbits there is a minimum impulse flight plan consisting of 
two impulses of short duration; see Fig. 5. The first impulse 
changes the vehicle’s orbit from the initial circle to the trans- 
fer ellipse which is tangent to the final desired orbit; the 
vehicle proceeds in free flight along the ellipse. The second 
impulse is applied at the point of tangency to the final orbit 
and the vehicle thereafter moves along the final circular orbit. 
Each impulse must be applied over a time interval that is 
very small compared to an orbit revolution. For missions 
originating from a low altitude above the Earth, Hohmann 
flight plans are limited to chemical propulsion. 

Irving and Blum (1) and Stuhlinger (4) have made analyses 
of transfers between circular orbits when using power limited 
low thrust propulsion systems. Examples of these are elec- 
trical or solar heater (3) propulsion. The thrust acceleration 
of these systems is much less than 1 g (32 fps), and the 
Hohmann flight plan with two short intense impulses cannot 
be achieved. These authors show that the shortest flight 
time is obtained when the thrust is applied continuously, 
producing a spiral orbit. However, they did not incorporate 
the possibility of storing some of the energy, and in this re- 
spect, their analyses were incomplete. It can be shown that 
the total thrust impulse required to make a transfer between 
two orbits depends upon the flight plan, the total impulse for 
the Hohmann flight plan being much less than that for the 
spiral orbit flight plan. With the use of energy storage flight 
plans can be obtained that have the same total impulse as 


An original version of this paper was presented at the ARS 
13th Annual Meeting, New York, Nov. 17-21, 1958, 


titled ‘‘Reduction of Flight Time and Propellant Requirements of 
Satellites with Electrie P 
cal Energy.” 

1 The work was sponsored in part by the U. 8. Air Force Office 
of Scientific Research, ARDC, under contract AF 49 (638)-61. 

2 Numbers in parentheses indicate References at end of paper. 


-ropulsion by the Use of Stored Electri- 


Use of Energy Storage in Low 
Thrust Spaceflight 


programming the thrust while still utilizing all the available energy. 


short thrust impulse at a point along the orbit. 
and used repeatedly. The performance with and without energy storage is compared for two types 
of missions, both in the gravitational field of the Earth: Round trips from small to larger circular 
orbits, and rendezvous maneuvers between orbits at low altitudes. 
with the extra weight of the storage batteries, flight plans can be achieved in some cases which re- 
duce both the flight time and the propellant requirements. 
that thrust be applied at only limited portions of the orbit. 
_ energy can be utilized, thereby substantially decreasing the flight time and the guidance problem. 
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is _ ‘This paper analyzes flight plans for spacecraft operating at thrust levels that are much smaller 
i¥ than the local gravitational force, e.g., electrical propulsion near the Earth. The thrust is limited 


For power limited flight, it has been suggested by Irving and 


Blum that the shortest flight times are obtained when operating at constant (thrust) acceleration. 
However, there was no consideration of the use of energy storage which permits more freedom for 


In particular, the electric 
over a portion of an orbit revolution and used to provide a 
The batteries would be recycled each revolution 


It is shown in case 1 that even 


For case 2 many flight plans require 
With energy storage all the available 


for the Hohmann flight plan. This can lead to a reduction 
in both the propellant consumption and flight time. 

This paper analyzes flight plans in which a part of the 
energy is stored for a portion of the orbit flight. This energy 
is then used more effectively when expended for short thrust 
impulses at desirable points along the orbit. The storage 
system is recycled during each orbit, thereby permitting ¢ 
sequence of thrust impulses. A comparison between pro- 
pulsion systems with and without energy storage is presented 
for two missions within the gravitational field of the Earth: 
1 Flight from an altitude of 100 miles to an orbit at a larger 
altitude, and subsequent return to the Earth; 2 simple in- 
terception and rendezvous missions between two satellite 
vehicles at low altitudes. 

Energy storage is also advantageous for missions at low al- 
titudes. In most low altitude maneuvers (for example, trans- 
fers from circular to elliptical orbits, or between orbits that 
are not in the same plane), optimum operation is achieved 
when impulses are applied only over limited portions of the 
orbits. The use of energy storage in low thrust propulsion 
systems permits the application of short impulses while still 
utiJizing all available energy. This decreases the flight time 
for the maneuver, and considerably simplifies the guidance 
and control problem. 

As a specific example, the present paper considers the use 
of storage in electrical propulsion systems, for which electric 
energy is stored in batteries. The advantage of storage is 
that better flight performance is obtained with short thrust 
impulses at limited regions along the orbits, while utilizing all 
the available energy. The disadvantage is the weight penalty 
of the storage batteries that must be used. The extent to 
which electric storage is useful for reducing the flight time and 
propellant consumption depends mainly upon the relative 
weights of the batteries and the rest of the propulsion system. 
The practical use of batteries was examined for missions to 
lunar orbits. It is shown that with nickel-cadmium cells, 
battery storage is of advantage when the specific weights of 
propulsion systems are about 20 lb/kw or more. 
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Missions to Large Circular Orbits 


In this section missions that require transfers between low 
and high altitude circular orbits are considered. The flight 
plans with and without the use of electric energy storage are 
evaluated, and a comparison is made of the propellant require- 
ments and the flight times. In order to use a minimum of 
calculations, two approximate flight plans are treated since 
the results can be expressed in simple analytical form: 1 A 
series of elliptical orbits produced by short thrust impulses; 
2 the spiral orbit resulting from continuous thrust. 

Some preparatory material is now presented with which 
this evaluation is made. First, satellite orbit nomenclature 
and kinematics are described. Next, the total mission im- 
pulse for a Hohmann transfer ellipse, and for elliptical and spi- 
ral orbit fight plans are determined. Then the formulas for the 
payload weight ratio are presented as a function of the flight 
time for given mission parameters. Finally, a comparison 
is made of the performance with and without energy storage 
using these preparatory results. 


Orbit Description and Kinematics 


We will assume that the only gravitational force felt by the 
satellite is the inverse square force acting from the center of 
the Earth. Perturbations due to the equatorial bulge, air 
drag and other celestial bodies are neglected. Thus, the orbit 
is an ellipse which is described by the variables shown in 


Fig. 1, where 
R = distance from Earth’s center to satellite 
R, = distance from Earths’ center to perigee point 
R, = distance from Earth’s center to apogee point 
6 = angle that radius vector makes with perigee line 


The equation of the ellipse is 
a(1 — e?) 


1 + ecos 0 


R= 
where the semimajor axis a = 3(R, + R,) and the eccentricity 
e = (R. — R,)/(Ra + R,). In terms of R, and R, 


wine (Ra + Rp) + (Ra — Rp) cos 0 [2] 


Also shown in Fig. 1 is the Earth with a mean radius Re. 
The altitude of the satellite above the surface of the Earth is 
given by h = (R — Rx). The apogee and perigee altitudes are 
ha = (Ra — Re) andh, = (R, — Re), respectively. 
The period of an orbit revolution 7 is given by the formula 
6) 
444 


V 


where go (32.1 ft/sec?) is the gravitational acceleration at the 
surface of the Earth. Fig. 2 presents the orbit period as a 
function of altitude h for two cases: Circular orbits, and el- 
liptical orbits with perigee altitude at 100 nautical miles and 
apogee at h. 

In order to determine the velocities of the satellite in a given 
orbit, it is convenient to start with the equation for the total 
energy which is constant in time. The total energy E is equal 
to the potential plus kinetic energies as shown in the formula 


[4] 


R 2 


Note that the total energy is zero when the satellite is at rest 
a very large distance from the Earth. Constant angular mo- 
mentum along the orbit gives the result that Rava = R,v, 
(6). Inserting this relation into the equation for the energy, 
we obtain the desired form for the velocity of the satellite at 


JANUARY 


SATELLITE 


Fig. 1 Notation used for the description of the elliptical orbit 
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Fig. 2 The period of orbit revolution for circular orbits at 
altitude 4 and for elliptical orbits with perigee at 100 nautical 
miles and apogee at altitude / 
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any point along the orbit 


1 
v= 2goRr? (; 


* 
Thus the velocities at apogee and perigee are 


) 
Ra + Rp 


| 


Va =e! [6a] 
(Ra + R,)Re Vint RR, 
For the special case of a circular orbit the velocity v, is 
R Rr + h 


Equation [6b] is plotted in Fig. 3 and denoted by v,.. The 
other two curves give the velocities at apogee va and perigee 
v> for an elliptical orbit with perigee at 100 miles and apogee 
at an altitude h as calculated from Equation [6a]. Also 
shown are three intervals Av,, Ave and Av,. These intervals 
have the following physical significance: Av, is the difference 
in velocity at 100 miles altitude between a 100-mile circular 
orbit and an elliptical orbit with perigee at 100 miles and an 
apogee altitude at h. Similarly, Av, is the difference in 
velocity at an altitude h between a circular orbit at an alti- 
tude h and an elliptical orbit with apogee at A and perigee at 
100 miles. Av, is the difference of velocity between two 


circular orbits at altitudes 100 miles and h. Fig. 4 shows the 
variation of the incremental velocities Av,, Ava, Av, and 
(Av, + Av.) with h as just described in Fig. 3. 
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ss 40,000 Total Impulse for Missions 

a 35,000 In order to obtain simple and analytic results, the total 
_ impulse requirements and the flight times will be determined 
ue 30.000 for two extreme flight plans: 1 Short intense thrust im- 
z pulses; 2. spiral orbits. In case 1 the thrust is assumed to be 
e; 25,000 much larger than the local gravitationa] force, and the latter 
- is neglected in determining the change in the orbit during the 
© 20,000 short impulse. Because of this assumption of rectilinear mo- 
= tion during the short impulse, the change in the vehicle’s 
> 15,000 velocity is equal to the integrated thrust acceleration. For 
w case 2 of the “tightly wound” spiral flight plan, the thrust is 
F 10,000 t smaller than the local gravitational force, and is applied con- 
- tinuously in the direction of motion. The orbit, if initially 
wt 5,000 circular, remains essentially a circle but of varying radius 
= ve during the mission. Note that in a circular orbit a satellite’s 


oan io? it io io speed remains constant due to the balance of the gravitational 
h-ALTITUDE IN NAUTICAL MILES and centrifugal forces. Thus, for the tightly wound spiral 

Fig. 3 0 chester veloc- where the orbit is essentially always circular, any change in the 
ities vp and va at perigee and apogee, respectively, in an elliptical vehicle’s speed is equal to the velocity change produced by 
the thrust chamber. Thus, for these two extreme cases, we 


orbit with perigee at 100 miles and apogee at an altitudeh 
will make use of the simple result that the velocity increment 


- produced by the integrated thrust acceleration during the transfer 
Pe i between two orbits is equal to the difference in speed of the 
ioe = vehicle in the two orbits. By assuming in case 2 that the 
 vehicle’s orbit remains essentially circular throughout the 


35,000 - + thrust period, the vehicle’s gravitational potential energy is 
o | 7 always twice the magnitude of the kinetic energy; see Equa- 
Y 30,000 ' tion [4]. Any change in the kinetic energy requires that the 
energy have twice as much change. Thus, it can be 
- 25,000 | eae shown that when applying positive thrust the vehicle’s velocity 
= Av, | decreases, and that this velocity change is equal in magnitude 

S 20,000 a 1 _ to the integrated thrust acceleration but of opposite sign. 

~ 15,000 Hohmann Flight Plan 

e The Hohmann transfer ellipse flight plan is presented here 
| because it represents the minimum total impulse plan. While 
2 ies! ‘it cannot be achieved with low thrust systems, it will serve as 
a a guide in determining other flight plans. The two impulses 
¥ required for a transfer from an initial circular orbit at 100 
ff miles toa large circular orbit are shown in Fig. 5. For very 
short intense impulses at perigee and apogee, the required 
h- ALTITUDE IN NAUTICAL MILES velocity changes of the vehicle are the Av, and Av, given in 
Fig. 4 Velocity differences as defined in Fig.3 Fig. 4. As shown in the previous paragraph these velocity 


changes are equal to the integrated thrust acceleration that 
the thrust chamber must supply for this mission. Since the 
Hohmann transfer ellipse requires the minimum total im- 
pulse, the sum (Av, + Av.) is the minimum velocity change 
due to thrust for a transfer between two circular orbits. 


Spiral Orbit Flight Plan 


In the spiral orbit flight plan the thrust is applied con- 
tinuously for the transfer between two circular orbits. As 
shown previously in the limiting case of a tightly wound 
spiral, the integrated thrust acceleration is equal to the 
difference in velocity of the two orbits. This velocity incre- 
ment is Av, for a mission from 100 miles altitude to a large 
circular orbit; see Fig. 4. Note that Av, is always larger 
than the sum (Av, + Av). Thus, other factors remaining 
unchanged, the spiral orbit flight plan always requires more 
propellant consumption than the Hohmann transfer ellipse. 
In fact, it is this difference that leads to the advantages of 
electrical energy storage for such missions. 

It is interesting to discuss the reason why the total impulse 


FINAL ORBIT requirements vary with the flight plan in making a transfer 
ysis between two given orbits. A vehicle moving along a given 
To tye orbit in space has a total energy that is independent of the 


may: ator - way in which the vehicle reached that particular orbit. Con- 
sider a mission between two orbits of energies and To 
ori effect a transfer between these two orbits would require the 
— irrespective of the flight plan 


ONE WAY: Av=(Avot Avg) 
ROUND TRIP: Av= (Av, + 2Av,) 
Fig. 5 Flight plan for the Hohmann transfer ellipse = (F-v)dt [7] 
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where F' and v are the exhaust thrust force and the vehicle 
velocity vectors, respectively. We will assume that the 
thrust is applied in the direction of motion, thereby producing 
the most work. The vector equation becomes scalar and the 
integral reduces to 


E, — E, = S Fedt [8] 
The usual definition of the average velocity 
S wat/S gat [9] 


where g is an arbitrary function. In this case S Fodt = 
oS Fat. Note that S Fat is the total thrust impulse for the 
mission and is equal to 


S = (E, [10] 


Since the energy difference is independent of the flight path, 
the lowest total impulse is obtained when the thrust is applied 
at the highest satellite velocity. In more general terms, the 
thrust should also be applied in the direction of motion at the 
point of lowest potential energy. The Hohmann transfer 
ellipse flight plan meets these requirements for the first im- 
pulse. There are several ways to reach the final circular 
orbit by using one or more impulses. The Hohmann transfer 
ellipse is best only when both the impulse and the flight time 
are considered. 


Elliptical Orbits With Low Thrust Systems 


Electrical propulsion systems are limited to low thrust ac- 
celeration. For example, to obtain an acceleration of 10~*go 
with specific impulse of 1000 sec, the total vehicle weight 
must be less than 5 lb for each kilowatt of directed power in the 
exhaust. This weight would be considered low for the power 
supply alone. Even for such a thrust operating continuously 
in the direction of motion, the shape of a low altitude orbit 
cannot be appreciably distorted. Starting with a circular 
orbit, a spiral orbit would result. 

A method for approximating the minimum total impulse 
flight plan with the use of a low thrust propulsion system will 
now be presented. It is based on the fact that with the 
provision for energy storage, all the electric energy produced 
can be applied in short impulses. Consider a mission from 
100 miles altitude to a high altitude orbit. Electric power 
produced at a constant rate (using solar or nuclear energy) 
would be stored in batteries over a portion of an orbit revolu- 
tion. In the vicinity of perigee, the thrust chamber would 
be turned on and the batteries discharged producing a short 
impulse.* The resulting flight path is shown in Fig. 6. With 
the short impulses at perigee, the apogee altitude would in- 
crease, but the perigee altitude would not change. This 
series of small thrust impulses at perigee requires the same 
total impulse as the single impulse Av,, used in the Hohmann 
transfer ellipse. When the apogee reaches the desired alti- 
tude, then impulses are applied at this apogee point in order 
to increase the perigee and thereby circularize the orbit. 

The advantage of the elliptical orbits flight plan (using 
storage) over the spiral orbit plan (continuous thrust) be- 
comes greater for a round trip from a low to a high altitude 
orbit and return. For the elliptical orbits flight plan part of 
the return trip can be assisted with atmospheric drag: The 
apogee pulses would make an elliptical orbit whose perigee 
lies inside the Earth’s atmosphere. The collapsing of this 
ellipse into a low altitude circular orbit can be accomplished 
with atmospheric drag on the vehicle. 


’For most electrical propulsion systems (are jet, magnetic 
accelerators, but not the ion rocket), the higher thrust level re- 
sults in a more efficient operation of the thrust chamber. 
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Fig. 6 Flight plan for a satellite using electrical propulsion with 
energy storage to go from a low altitude circular orbit at / to a 


larger circular orbit at ¥. Short impulses are initially applied 

at the point / producing a series of expanding ellipses 1, 2, 3 and 

4. When the apogee of the expanding ellipses reaches the final 

circular orbit then pulses are applied at this apogee point A 

thereby increasing the perigee altitude and circularizing the 
ellipses 5, 6 and 7 


=2AVs 


fe) 
8 


AVeLuiptic Avg) 


je) 


10,000 +-——— 


1d 10 10 
h- ALTITUDE IN NAUTICAL MILES 


AV -VELOCITY CHANGE IN FT/SEC 


Fig. 7 Velocity increments Av required to go from 100 miles 
altitude to a circular orbit at altitude h and return to Earth. 
Adspirat ANd Aveiliptic are the velocity increments for the spiral orbit 
expanding ellipses flight plan, respectively 


Fig. 7 shows the velocity increment Av required for a round 
trip from a 100-mile circular orbit to a larger circular orbit of 
altitude A and return. The curve denoted by Aveniptic is for 
the minimum impulse flight plan shown in Fig. 6 with at- 
mospheriec drag circularizing the return ellipse. The other 
curve Adspirai is for the spiral orbit flight; atmospheric drag 
cannot be used. The formulas for these velocity increments 
are 


= 


Aveniptic = (Av, + 2Ava) 


Abspiral = [11] 
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Practical Propulsion Systems Using Battery Storage 


The use of energy storage requires the incorporation of elec- 
tric batteries in the propulsion system. Since these batteries 
have to be recycled many times the battery materials and the 
mode of operation are limited. Nickel cadmium and silver 
cadmium cells appear to be the most practical at the present 
time. These batteries have an operating weight of about 20 
to 25 lb/kw-hr storage capacity. In order to maintain reliable 
operation over many recycles, the charging time must be at 
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CHARGE 
POWER DIRECT 
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30 MINUTE 
BATTERY DISCHARGE 


Fig. 8 Thrust program when using electrical energy storage for 
expanding the orbit 
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Fig.9 Instantaneous apogee h. and perigee h» altitudes variation 
with time for the particular example given in text. The satellite 


makes many elliptical revolutions similar to the orbits shown in 
Fig. 6 


least 1 hr and the discharge should be at least 3 hr. The 
newly developed fuel cells show promise for a lighter weight 
with higher storage capacity. In the present discussion we 
have considered only the commercially available batteries. 
A practical compromise has to be made between the most 
desirable storage capacity and the weight penalty of the 
batteries. For low altitude circular orbits, the orbit period is 
of the order of 1.5 hr. Since the batteries are recycled each 
revolution, the storage time need not be larger than the period 
of revolution. For larger orbits the orbit period increases 
with increasing apogee as shown in Fig. 2. However, the 
period is still only 2 hr even when the apogee is at 2000 
miles altitude. When expanding into a highly elliptical orbit, 
the major part of the thrust is required in the initial expansion 


of the apogee from 100 miles to several thousand miles alti- 


tude, because the gravitational forces decrease with the 
square of the distance from the Earth’s center. 

A practical flight plan is shown in Fig. 8 where a 2-hr 
storage capacity is included in the electrical propulsion sys- 
tem. When the period of revolution is less than 2.5 hr, then 
the left figure applies. The solid line shows the free flight 
part of the orbit during which the electric energy production 
goes to the batteries for storage. The batteries are discharged 
into the thrust chamber for 30 min around perigee. These 
charging and discharging rates are consistent with long-life 
operation of the batteries. For periods greater than 2.5 hr the 
figure to the right applies; for 2 hr around apogee the electric 
energy is stored in the batteries, and for 30 min around perigee 
the batteries are discharged. During the rest of the orbit the 
battery system is not in use and the power goes directly to the 
thrust chamber. In this way all the energy produced is used 
for thrust. 

We have carried out exact calculations (using a high speed 
computer) in order to determine the flight path when trans- 
ferring from a 100-mile circular orbit to a lunar orbit using the 
thrust program similar to that described in Fig. 8 with 1, 1.5 
and 2-hr battery storage capacity. Some of the results are 
presented in Fig. 9 for a 2-hr storage capacity case. The 
curves show the apogee and perigee altitudes as a function of 
time; the vehicle actually makes many elliptical revolutions 
similar to the orbits shown in Fig. 6. For this example, the 
total thrust requirement was within 10 per cent of that for the 
Hohmann transfer ellipse. Thus, little is to be gained in 
more than 2-hr storage capacity. When using 1-hr storage 
capacity, the thrust requirement was about 50 per cent 
greater. Different conclusions would be obtained when con- 
sidering a transfer to another altitude. Obviously, the mis- 
sion requirements and the specific weight of the power supply 
must be known before the amount of storage can be optimized 
and the advantages of storage can be evaluated. 


Optimizing Payload Weight Ratio 


In order to make a comparison of the spiral and elliptical 
orbit flight plans for specific missions, the payload weight 
We will use the 


ratio and the flight time must be calculated. 
following notation 


a = thrust acceleration 
exhaust velocity 


P = thrust power in the exhaust : 

k = structure factor (propellant container is assumed pro- 
portional to the propellant weight ) 

m = vehicle weight 

mp = battery weight 

m; = initial vehicle weight 

= payload weight 

Mp, = initial propellant weight 

Mw» = power supply weight 

T = mission flight time 

ve = V27(1 + k)/a = characteristic velocity 
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My/P = 


mp/P = 


specific weight of propulsion system 
specific weight of batteries 


A given mission will have parameters that are fixed while 
others can be optimized. For example, the fixed quantities 
may be flight time 7; velocity change due to thrust Av; 
structure factor k, and specific weight a. The expression for 
the payload weight ratio (mz/m;) will be derived assuming 
that the exhaust power P, averaged over an orbit revolution, 
is constant. Conservation of mass gives the relation 


(1 + + my + [12] 


m = 


when the subscripts 7, p, w and L are defined in the foregoing. 
The use of the structure factor k assumes that the propellant 
container is proportional to the propellant. The propellant 
required for the mission is 


= fi —d = — [13] 
0 Qa 0 
+. 
c = exhaust velocity oe 
@ = specific weight of the propulsion system 
= total flight time 
Integrating the rocket reaction equation a = me/m 
mp = — So [14] 


where a and ¢c are the instantaneous acceleration and exhaust 

velocity. These equations combine to give the payload weight 

ratio 

m; of (dt/c?) 
0 

No attempt has been made to maximize the payload weight 
ratio for all possible flight paths between two circular orbits. 
However, the results presented are for two important cases: 
1 Constant acceleration, and 2 constant exhaust velocity. 
It has been shown (1) that a larger payload weight ratio is 
obtained for the constant acceleration case. However, the 
difference between the two cases is negligible when the pay- 
load weight ratio is over 20 per cent as shown in Fig. 10. Case 
2 with constant exhaust velocity is applicable to systems like 
the nuclear rocket or solar heater where it is desirable to 
operate at the maximum specific impulse obtainable. 

Cask 1: The acceleration a is kept constant. The payload 
weight ratio is obtained by first eliminating c in the equations 
me = am and mc? = 2m./a, and then integrating. One ob- 
tains the relation 


m,; a( Av)? 


[16] 


Mp 


where Av is the velocity increment produced by thrust ac- 
celeration. Substituting for m, in Equation [8] the payload 
weight ratio becomes 


(27 /a(Av)?] + 1 


my 


My 
= [17] 
mM; 


where (m,,/m;) can take a range of values. mz ,/m, is a maxi- 


mum when 


Ve Ve 


where v, is called the characteristic velocity and is defined 


[18] 


270 + k) 


[19] 


= 
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‘ adt = velocity increment due to thrust acceleration 
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Fig. 10 The maximum payload weight ratio as a function of 

Av/c. The structure factor k was fixed at 0 and 0.1. The solid 

curves are for constant acceleration and the dashed curves are 
for constant exhaust velocity flight plans 


The maximum payload weight ratios becomes 


\2 
mi Ue 


which is only a function of Av/v, and k. 


Case 2: The exhaust velocity ¢ is kept constant. Then 
Equation [15] becomes 
2 
— (1 — [21] 
Ue 


~The exhaust velocity c is the only variable, as Av, k, a and v- 
fixed quantities. 


The maximum value of mz,/m; is ob- 
tained when the following relation is fulfilled 


Note that (Av /c) is in implicit function of (v./Av) only. Solv- 
ing Equation [22] for (Av/c), and inserting it in Equation 
[21] leads to an expression for the payload weight ratio 
which is only a function of (v./Av) and k. Fig. 10 shows a 
plot of the maximum payload weight vs. (Av/v.); k was fixed 
at 0 and 0.1. Equation [20] is used for the constant accelera- 
tion case, and Equations [21 and 22] combined for the case 
of constant exhaust velocity. It should be noted that the 
curves practically overlap when (mz,/m;) is greater than 0.2. 
Thus for practical flight problems there is no essential ad- 
vantage of the constant acceleration over the constant thrust 
flight plan; in fact, operation of thrust devices over too wide 
a specific impulse range may require inefficient operation. 

For the case of constant exhaust velocity the maximum 
payload can be expressed in terms of k and Av/c by substitut- 
ing for v, from Equation [22] into Equation [21] 


1 — [1 + 


mL 


mi 


Equation [23] is plotted in Fig. 11. It is not possible to obtain 
a flight plan that could be represented by a point in the 
region to the left of these curves. Thus, such curves can be 
used to determine the optimum exhaust velocity ¢ for each 
mission; i.e., for a given velocity increment Av and desired 
payload weight ratio. 


Comparison of Flight Time and Payload 


The framework necessary for comparing the performance of 
electrical propulsion systems with and without electrical 
storage has thus been prepared. We will use two of the 
several ways in which this comparison can be made: 1 The 
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Fig. 11 The maximum payload weight ratio plot as a iden : 


of (c/Av) for constant exhaust velocity c. Av is the velocity incre- 
ment that the thrust acceleration must supply for the mission. 
The structure factor & was fixed at 0 and 0.1 Papi 
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Fig. 12 The ratio of the flight times for the same payload weight 

ratio as a function h for a trip from 100-mile altitude to a circular 

orbit at A and return to Earth. a and az are the specific weight 
of the propulsion system and the storage batteries 
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Fig. 13 The payload weight ratio as a function of flight time for 

a round trip from 100 miles altitude to the 24-hr orbit (22,000 

miles altitude). Two values for the propulsion system specific 

weight are shown: a = 10 lb/kw and a = 60 lb/kw. The 
battery specific weight az = 30 lb/kw 
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flight times are compared for a fixed payload weight ratio; 


2 the payload weight ratio is compared for the same flight 
times. 


In the previous section it was shown that the payload weight 
ratio is only a function of (Av/v,) and k. Thus, for a fixed 
payload weight ratio mz/m; and structure factor k, (Av/v,) is 
constant independent of the flight plan. This important re- 
sult states that the flight time 7’ is given by the relation 


Ade 


for all flight plans. Thus, the ratio of the elliptical and spiral 
orbit flight times is 


T ciiiptical (Geniptieat) (Aventiptical)® (25] 


The structure factor is assumed to be the same for both cases. 
The Av required for a round trip mission to a circular orbit 
and return to Earth is shown in Fig. 7. Avenipticeaa = (Avy + 
2Ava) and Adspira: = 2Av,. For the elliptical orbit, the extra 
weight of the batteries, say mg, must be added. The weight 
of the whole propulsion system becomes m, + msg. The 
corresponding specific weight Qenipticar = (m,/P) + = 
a + az where az is the specific weight of the batteries alone. 
For the spiral orbit case there is no storage sO Qspira1 = @ 


assuming the same propulsion system in both cases. Sub- 
stituting these values for Av and a in Equation [24] 
Teiiptical _ & + Qp + 2Ara\? [26] 
piral Qa 2A, 


The specific weight of the batteries ag for a practical design 
is estimated to be from 20 to 40 lb/kw depending on the 
mission. The power supply specific weight will vary between 
10 and 60 Jb/kw depending on the size of the power supply. 
Thus the ratio (@z/a) will range from 0.5 to 3. In Fig. 12 is 
shown the variation of the flight time ratio given by Equation 
[26] for round trip missions from a 100-mile circular orbit to 
larger circular orbits. The curves are for (ag/a) equal to 3 
and 0.5. Note that when ag/a = 0.5, corresponding to a =~ 
60 lb/kw, the flight time is always less for the elliptical orbits. 
For az/a = 3 (a = 10 lb/kw), elliptical orbits become more 
favorable only above 40,000 miles altitude. 
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Fig. 14 The payload weight ratio as a function of flight time for 

a round trip from 100 miles altitude to a lunar orbit. Two values 

for the propulsion system specific weight are shown: a = 10 

Ib/kw and a = 60 lb/kw. The battery specific weight ag = 
30 lb/kw 
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The payload weight ratio as a function of the flight time is 
shown in Figs. 13 and 14 for two missions: A round trip from 
a 100-mile circular orbit to the 24-hr orbit (22,000 miles alti- 
tude) and lunar orbit. The curves were obtained with the 
use of Equation [20]. The (Av)’s are for the spiral and ellip- 
tical flight plan from Fig. 7. Note that the payload weight 
can always be exchanged for flight time. ae v7 


Limitations on This Simplified Analysis 


The analysis presented above to evaluate the usefulness of 
electrical storage was simplified considerably. The justifica- 
iion was that simple analytical results could be obtained 
which are in fair agreement with the practical flight situation. 
However, this analysis for the spiral orbit flight plan becomes 
very poor when the local gravitational force is equal to ov less 
than the thrust force; the flight changes from a spiral to a 
straight line (7,8). This occurs at an altitude of approxi- 
mately 80,000 miles for a thrust acceleration of 10~* go and 
at the altitude of the lunar orbit 10~4 go. However, even for 
a 10~* go acceleration, the Av required to a lunar orbit is still 
the Av, shown in Fig. 4, because the change from tightly 
wound spiral orbit is compensated by the extra thrust re- 
quired to return to a circular orbit. This compensation does 
not hold for the return trip, and for 10~* go the Av require- 
ment becomes less than the Av, of Fig. 4 for altitudes above 
100,000 miles. For 10~* go our analysis cannot be used for 
altitudes above about 200,000 miles. 

In a practical situation the battery storage capacity is 
limited, and it is not possible to achieve the ideal flight plan 
with impulses at perigee and apogee, as shown in Fig. 6. 
The thrust requirements and the flight time will thus be 
larger than that given in our simplified analysis. It was shown 
that in a practical situation about 10 per cent more flight 
time is needed for trips to a lunar orbit. This would be less 
for low altitude orbits. On the other hand, the spiral orbit 
flight times are over a factor of 2 longer than the ideal ellipti- 
cal case. Thus, our analysis is a valid approximation. 


Summary 


The usefulness of electrical energy storage for effecting a 
transfer between two circular orbits has been evaluated. 
The battery weight penalty required for the electric storage is 
compensated by the smaller thrust requirements of the ellipti- 
cal flight plan. In evaluating a particular mission an impor- 
tant parameter is the relative weights of the batteries and the 
rest of the propulsion system. For batteries that are pres- 
ently available the electrical storage becomes practical for 
some missions when the specific weight of the electrical pro- 
pulsion system is greater than 20 lb/kw. Another ad- 
vantage of the use of the elliptical orbit flight plan is available 
for the initial portion of the trip where the perigee altitude 
remains low. If, for some reason, it becomes necessary to in- 
terrupt the trip, a return to a lower orbit or to the Earth can 
be accomplished quickly. An impulse is made near apogee 
which lowers the perigee altitude into the atmosphere. The 
orbit can be circularized quickly with the use of a suitable at- 
mospheric drag system. 


Interception and Rendezvous Problems 


In this section we will consider missions between orbits in 
which the total energy of the vehicle does not change appre- 
ciably, that is, transfers between orbits at essentially the 
same altitude above the Earth. These flight plans arise in 
interception and rendezvous missions. It is convenient to 
consider separately the motion along three coordinates: 1 A 
transfer between two orbits that are in the same plane, 2 a 
mission requiring a change in latitude 7 of the orbit plane, and 
3 a change in the angle Q, i.e., a rotation of the line of nodes 
of the orbit plane. 
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After the impulse, the change in period becomes 


ORBIT OF A AFTER 
Ist IMPULSE 


1% 


Fig. 15 The flight plan to make a rendezvous between two 
satellites A and B that are moving on the same orbit but at 
different positions 


Transfer Between Co-Planar Orbits 


The problem of determining the best flight plan for making 
a rendezvous (or a transfer) between two arbitrary orbits at 
essentially the same altitude is more involved than the Hoh- 
mann problem of transfer between two circular orbits. We 
have only considered the simpler problems of: 1 A ren- 
dezvous between two vehicles moving along the same ellipti- 
cal orbit but at different points on the ellipse; 2 a transfer 
from a circle to an ellipse. The more general problem of 
transferring between two arbitrary ellipses has been con- 
sidered by Lawden (5). 

Fig. 15 shows the flight plan for the rendezvous between 
two vehicles A and B that are in the same elliptical orbit but 
separate by a fraction N of the orbit period. Suppose that 
vehicle A can maneuver and wants to make contact with 
vehicle B. This can be accomplished with two impulses. 
The first impulse is applied tangentially to the orbit with the 
exhaust thrust pointed in the direction of B, resulting in a re- 
duction in the flight energy of A. The orbit of A now changes 
to a smaller ellipse remaining tangent to the original ellipse at 
the point of impulse, at A. The period of revolution 7 for the 
resulting elliptical orbit is smaller than the original period. 
Thus, as time progresses for each revolution of satellite A, B 
will come closer to the tangent position of A. If the change in 
period Ar, of satellite A, is 1/n (where n is an integer) of the 
original separation time 7N, then after n revolutions of A, 
satellite B will recede to the tangent position of A. When A 
and B become coincident, then A makes a reverse impulse 
which gives it the same velocity as A and the rendezvous is 
completed. 

We will now calculate the propellant requirement for such 
a mission. The velocity of the vehicle at point A is 


vg? = 2gokz? EB = | [27] 
‘Ss 

where a is the semimajor axis. When an impulse producing 
the velocity change Av is applied at A, the semimajor axis a 
changes by Aa, but Ry, remains constant. The relation be- 
tween Av and Aa is 


ORBIT OF SATELLITES A AND B 
BEFORE Ist IMPULSE 
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Inserting Equation [30] into Equation [28] yields Av requirements can be very small by making Ar small and 
waiting a long time. This very small Av is possible because 
ey goRe? Ar 31 initially the satellites are in the same orbit. 
— a. a. [31] Equation [32] indicates that for a given ellipse the impulse 


requirement is minimized when v, is a maximum, since all 

If the satellites become coincident after n orbit revolutions, other quantities remain unchanged. This indicates that 
the time 7 necessary for the mission is 7 = n(r — Ar), and satellite A should wait until the perigee is reached and then 

the required period change Ar is Ar = Nr/n, where N is the apply the impulse. This is the most efficient way of apply- 
fraction of the period that the satellites are initially separated. ing thrust. ; : . ; 
Equation [31] becomes When using electrical propulsion without energy storage it 

is possible to accomplish the rendezvous if satellite A uses a 


stant retarding thrust for half the flight time and ac- 
goRr?t — Ar constan g 

Av = Wr [32] celerating thrust the latter half. If the acceleration is small, 
=| - ee . "Ae! the orbit will spiral in and then out although the orbit at any 


time will have the shape of the original ellipse. For this spiral 
orbit case the analysis has been made only for an initially 
circular orbit and the necessary impulse is twice that required 
for the previous case using two short impulses. 


The total impulse for the mission will be twice the one given 
by this equation, since a similar impulse has to be applied to 
return satellite A to its original orbit. It should be noted 
that the range of 7 in Equation [32] is not continuous, but 
oa | only be integral multiples of the period (r — Ar). The 


Transfer Between Circular and Elliptical Orbits 


The more complicated problem of making a rendezvous be- 
tween a circle and an ellipse can, in general, be separated into 
parts that have already been treated: 1 Making the first 
impulse of the Hohmann transfer ellipse; 2 the shifting of 
angular position along an orbit. This is not the only possible 
flight plan and there was no attempt to determine the opti- 
mum thrust plan. 


pa! 
——- The use of electric storage has a very large advantage for 


missions requiring a transfer from a circular to an elliptical 
Oey wr tf orbit. Since impulses are required at only limited portions 
= of the orbit, the flight time would be reduced with the use of 

ea storage. Flight plans as described in Figs. 6 and 8 would be 
used, 

Ji 

Baie > Missions Requiring a Rotation of the Orbit Plane 

4 i In order to rotate the orbit plane thrust must be applied 


perpendicular to the orbit plane. An instantaneous velocity 
- inerement Av applied to the satellite normal to the plane of 
the orbit will rotate the orbit plane about the radius vector to 
that point of an angle AX where 


v9 V goRz%a(1 — e?) 


— vg is the azimuthal component of the satellite velocity. The 

; : relation on the right side expresses vg in terms of the usual 
Fig. 16 Rotation of the orbit plane by an impulse Av applied orbit variables. This rotation is indicated in Fig. 16. An 
normal to the plane impulse applied at another point will cause a rotation about 


that point. 


M a a oan A small angular rotation AX can be resolved into rotations 
ean about the coordinates Q and 7, where Q is the angle between 
the vernal equinox and the line of nodes in the equatorial 
2 


a ; plane, and 7 is the inclination of the orbit with respect to the 


ORBIT PLANE ROTATION AX = AQ + Ai 


where 


i i 
sin 7 V goR —e?)sint 
THRUST 


IMPULSE Ai = Ad cos u = Av _—_—s«O[86 


= 
_ The angle between the radius vector to the satellite and the 
a __ Jine of nodes in the equatorial plane is designated by u. To 


= a produce the largest change in latitude angle 7, without any 
ree - change in ©, the impulse Av should be applied when passing 
through the equatorial plane, i.e., for uw = 0 and 180 deg. On 


THRUST 
IMPULSE 


Fig. 17 Thrust program when using electrical energy storage for the other hand, an impulse applied when u = +90 deg would 
rotation of the orbit produce the largest change in Q with no change in 7. 
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The best control when changing the direction of the orbit 
plane is accomplished with short impulses at desired points 
in the orbit. Such maneuvers which require pulsing are more 
suited to electrical propulsion having electric storage, since 
this permits the complete utilization of the available energy. 
The flight plan is shown in Fig. 17. The batteries would be 
charged during the free flight periods and discharged, produc- 
ing relatively large impulses, at the desired points along the 
orbit. Since the thrust pulses can be used with equal efficiency 
on opposite sites of the orbit, the storage capacity need be only 
about } of the period. This would correspond to less than a 
3-hr storage capacity for low altitude orbits. Thus, the re- 
quired specific weight of the batteries would be in the range 
from 10 to 15 lb/kw. It can be shown that there is an ad- 
vantage to the use of battery storage, even with this weight 
penalty, when the specific weight of the propulsion system 
(excluding the batteries) is greater than 15 lb/kw. The flight 
time and propellant requirements would be less. nits: 
Summary 


We have compared the performance of electrical propulsion 
systems with and without energy storage for several missions 
in the gravitational field of the Earth. The advantage of 
electric storage depends upon the relative weights of the 
electrical propulsion system and of the batteries, and upon the 
relative decrease in the impulse requirement for the elliptical 


orbit compared to the spiral orbit flight plans. 

Energy storage is advantageous for rendezvous maneuvers 
at low altitudes. For most low altitude missions (for example, 
transfers from circular to elliptical orbits, or between orbits 
that are not in the same plane), optimum operation is achieved 
when impulses are applied only over limited portions of the 
orbit. The use of energy storage in low thrust propulsion 
systems permits the application of short impulses while still 
utilizing all available energy. This decreases the flight time 
for the maneuver, and considerably simplifies the guidance 
and contro) problem. 
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Composition of Reaction Mixture 
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_ A simple procedure is developed which permits the pre- 
diction of near-frozen or near-equilibrium flow in rocket 
nozzles. It is found that the reaction 2 NO = N; + O2 
in aniline-nitric acid rocket motors is nearly frozen at the 
chamber exit, in agreement with the earlier finding of 
The recent data on the supersonic flow of N:O, 
The predicted com- 


Penner. 
in de Laval nozzles are examined. 
position agrees satisfactorily with the experimental value 
at the nozzle exit. 
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C' YNSIDER the reaction A — B, the kinetics of which is 
governed by the equation 


dx/dt = k(a — x) 


where 
k = rate constant 
a = initial composition 
x = amount of A decomposed 


rate constant is related to temperature by the relation 


The 
k = 
where 
E = energy of activation _ 
T = temperature 
C = certain constant 
R = gas constant 


The above equation can be integrated under isothermal con- 
ditions, and the composition of A can be predicted at any 
instant provided the rate constant is known. However, in 
practice, we come across situations where the temperature 
varies with time, and hence the above equation cannot be 
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integrated. Problems of this nature are encountered while 
estimating the composition of the reaction mixture across 
the rocket nozzle where the temperature suffers a sharp drop 
of several thousand degrees. The estimation of composition 
is essential for computing the correct exhaust velocity of 
— since the velocity v at the nozzle exit is given by 
(1/2)Mv H, — He {1] 


enthalpy of the fluid mixture at constant pressure when 
one-dimensional adiabatic flow is assumed 


The subscripts c and e denote chamber and nozzle exit, respec- 
tively. 

The composition of the reaction mixture at the exit is 
usually determined by assuming the freezing of chemical reac- 
tion at the chamber exit or by assuming the establishment 
of immediate chemical equilibrium at the nozzle exit (1).' 
For this purpose, useful criteria have been developed by 
Penner, but the greatest disadvantage is that in actual appli- 
cation the frozen flow criterion and equilibrium flow criterion 
have to be applied separately. 

The aim of the present paper is to evolve a more straight- 
forward procedure for deciding the nature of the flow in con- 
ventional rockets. 


Near-Equilibrium and Frozen Flow in Rocket Nozzles 


We shall consider, for the sake of simplicity, reactions of 

the type 

ky 

= Y+2Z2 

kb 
ky is the rate of forward reaction, and k, is the rate of back 
reaction. Reactions of the above type occur in aniline-nitric 
acid and hydrogen-fluorine rocket motors. We assume that 
at the chamber exit the number of molecules of the species 
X, Y and Z are a, b and ¢, respectively. If V is the volume 
and z is the number of molecules of species X decomposed, the 
rate of reaction is given by 


d 


Since ky and k, vary with temperature, the above equation 
cannot be integrated. However, some approximate con- 
clusions can be reached if it is assumed that the cooling rate 
is constant. Before doing that we shall discuss the general 
case. 

Since z is a function of temperature 7, time ¢ and pressure 
P, we can write by the methods of partial differentiation 


dt \at oP/r. aT} dt 


Here (—dT7'/dt) is the cooling rate which would depend on the 
dimensions of the nozzle. It can be seen that the net reaction 
rate can be determined if the quantities on the right-hand 
side of Equation [3] can be evaluated. This is done below. 


At equilibrium 
(a — 2) 
where 2, is the equilibrium value of zx. Using this relation we 


can write Equation [2] in the following form 


at ~ (b+ + 2/2) 


[5] 


3 Numbers in parentheses indicate References at end of paper. 
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dx dA dB 
( aT 


On integration, at constant temperature and pressure, this 
equation yields 
L 
Vi 4b + (c + 2/2) 


1 


2a%b + c) + Sabe + (a? — Abe)te 
(a? — 4bc) — 22a +b +c) 
3 


When ¢ = 0, x = 0, and hence Equation (6] reduces to 


1 Ieyt 1 (A + Ba)a, 
4b+2./2)(c+2-/2) A(t —2) 


from which we obtain 


— 1) 


where 
ky-(A + Bae) 
Differentiating both the right and left sides of Equation [8] 
with respect to 7, keeping time and pressure constant, and 


~ 4(b + + 
eliminating ¢ with the help of Equation [7] we have (2) 


— 
A(A + Ba.) 
(A +B 2) (a, 2) 
Be) 
(A + Ba)z. 


[9] 


Using the relation 


we have Sa. 
_ 2a — %) (b+ 2/2) (c+ 2/2) AH iy 


+6 +c) + (ab + ac + 4bc) RT? 


AH = heat of reaction a _ 


K = equilibrium constant 


From Equations [9 and 10], the value of (dx/dT);,p can be 
readily evaluated. In a similar way, one can evaluate 
(dx/dP) 7, and hence the net reaction rate can be computed. 

We can learn whether near-equilibrium or frozen flow is 
maintained by using the condition dx/dt = 0 at the exit, 
which would be true for both the cases. Thus, for both cases, 


we shall have 
dx aT 


dt dP) aT 


The resulting equation obtained after proper substitution can 
be solved for xz. If x is found to be zero, freezing of chemical 
reaction would take place. On the other hand, when z is 
equal to the equilibrium value, one can safely conclude that 
near-equilibrium flow occurs. 

As an illustration, we take the case of aniline-nitric acid 
rocket motor where the basic reaction involved is 


2NO=N:+ 


For the sake of simplicity, we shall neglect the intermediate 


(Continued on page 113) 
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